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Abstract. The line geometric model of 3-D projective geometry has the 
nice property that the Lie algebra sl( 4) of 3-D projective transformations 
is isomorphic to the bivector algebra of Cl( 3,3), and line geometry is 
closely related to the classical screw theory for 3-D rigid-body motions. 
The canonical homomorphism from SL( 4) to Spin( 3, 3) is not satisfy¬ 
ing because it is not surjective, and the projective transformations of 
negative determinant do not induce orthogonal transformations in the 
PKicker coordinate space of lines. 

This paper presents our contributions in developing a rigorous and conve¬ 
nient algebraic framework for the study of 3-D projective geometry with 
Clifford algebra. To overcome the unsatisfying defects of the Pliicker cor¬ 
respondence, we propose a group Pin 3p ( 3, 3) with Pin( 3, 3) as its normal 
subgroup, to quadruple-cover the group of 3-D projective transforma¬ 
tions and polarities. We construct spinors in factored form that generate 
3-D reflections and rigid-body motions, and extend screw algebra from 
the Lie algebra of rigid-body motions to other 6-D Lie subalgebras of 
si(4), and construct the corresponding cross products and virtual works. 

Key words: Projective Geometry; Line Geometry; Screw Theory; Pliicker 
model; Geometric Algebra. 


1 Introduction 

The study of the geometry of lines in space was invented by Pliicker with his 
introduction of the now so-called Pliicker coordinates of lines. It became an 
active research topic with the establishment of screw theory by Balls [1], where 
the 6-D Pliicker coordinates of a line are decomposed into a pair of 3-D vectors, 
called the screw form of the line, and the inner product and cross product of 
vector algebra are extended to screw forms. 

A pair of force and torque, called a wrench , are naturally represented by a 
pair of 3-D vectors, and are geometrically interpreted as a line in space along 
which the force acts, together with a line at infinity about which the torque 
acts. On the other hand, a pair of infinitesimal rotation and translation, called 
an infinitesimal screw motion or rigid-body motion or twist, are represented by 
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the rotation axis and the translation vector, and are again naturally represented 
by a pair of 3-D vectors. Geometrically the translation is a special “rotation” 
about an axis that is at infinity, so the translation vector represents a line at 
infinity. For a wrench (f, q) and an infinitesimal screw motion (v, t), where f is 
the force direction multiplied with the magnitude of force, q is the composed 
torque, v is the rotation axis direction multiplied with the angle of rotation, and 
t is the moment of the screw motion, the virtual work of the wrench along the 
infinitesimal screw motion is the “crossed” inner product 

(;)•(:) ! = f -* + —■ m 

The inner product (1.1) gives the 6-D space of wrenches a signature R 3,3 , 
where a pure force has zero inner product with itself, called a null vector. A 
positive (or negative) vector of R 3,3 is one having positive (or negative) inner 
product with itself. A positive vector is interpreted as a pure force together 
with an extraneous torque so that the pair follow the right-hand rule, while 
for the negative vector, the force and torque follow the left hand rule. The 
group SL( 4) which acts in the 4-space of homogeneous coordinates of points, 
can be lifted to a group action in the 6-D space of wrenches by acting upon 
the Pliicker coordinates of the lines representing the wrenches. The image of the 
lift is SOo(3, 3), the connected component of 50(3,3) containing the identity 
[3]. As PSL( 4) = 5L(4)/Z 2 is the group of orientation-preserving projective 
transformations, the crossed inner product provides a 6-D orthogonal geometric 
model of wrenches to study 3-D projective geometry of points. 

The same inner product (1.1) also gives the 6-D space of twists the same 
signature R 3,3 . The interpretation of a null vector of R 3,3 in the setting of twists, 
is that it represents an infinitesimal pure rotation or pure translation. A positive 
vector represents an infinitesimal screw motion where the translation along the 
screw axis follow the right-hand rule with the orientation of the rotation, while 
a negative vector represents a left-handed infinitesimal screw motion. The lift 
of group SL( 4) to 5C*o(3,3) then makes the space of infinitesimal rigid-body 
motions a 6-D orthogonal geometric model to study the orientation-preserving 
projective geometry of points. 

The bold-faced words clearly reveal a conflict. The group of rigid-body mo¬ 
tions is 6-D, while the group SL( 4) is 15-D; the former is much smaller. What 
sense does it make to investigate projective transformations via rigid-body mo¬ 
tions? Furthermore, the inner product (1.1) is between the space of wrenches 
and the space of twists, indicating that the two spaces need to be identified, 
yet they have to be different spaces by nature. Understanding (1.1) as a pairing 
between a linear space and its dual space does not make much difference, as the 
same inner product exists in either space. 

Line geometry and screw theory are closely related to each other. In history, 
the screw forms were first used by Clifford in the name of biquaternions, also 
known as dual quaternions, in describing 3-D Euclidean transformations. Later 
on, Balls [1], Study [35], Blaschke[4] established screw theory and developed dual 
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vector algebra out of Clifford’s dual quaternions, also known as screw algebra. 
Nowadays line geometry together with screw theory have important applications 
in mechanism analysis, robotics, computer vision and computational geometry 
[6], [2], [33], [10], [29], [32], [9], 

For two vectors (xi,yi) T £ R 3,3 and (x 2 ,y 2 ) T £ R 3 ’ 3 , where Xj,yj £ R 3 , 
their cross product , also called dual vector product, is defined as follows: 

( xi W X2 W xi ; x2 V (i.2) 

\yi ) \y2 J V x i x + yi x x 2 y 

This product is covariant under the subgroup of SOq( 3, 3) that is the lift of the 
group of rigid-body motions SE( 3), but not so under the whole group SOo(3, 3). 
In other words, it is not a valid operator in 3-D projective geometry; it is valid 
only for Euclidean geometry. 

In dual vector algebra, the dual inner product of two vectors of R 3 is defined 
to be a dual number. A dual number is of the form A + e/i where A, p, £ R, e 2 = 0 
and e commutes with everything. This numbers system is a ring instead of a 
field, and the corresponding polynomials and modules are drastically different 
from the usual ones. The dual inner product is invariant under the lift of SE{ 3) 
to SO 0 ( 3,3), so it is suitable for Euclidean geometry only. 

In [25], it was pointed out that dual vector algebra and dual quaternions can 
be realized in the conformal geometric algebra Cl( 4,1), and can be extended to 
arbitrary dimensions. The Euclidean geometric parts of the wrench model and 
the twist model have no conflict, and their identification is natural. The twist 
model should not have anything beyond Euclidean geometry, otherwise it would 
be absurd. The wrench model, or more generally the model of lines in space, 
deserves further attention. 

For invariant computing in projective geometry, the traditional algebraic tool 
is Grassmann-Cayley algebra and bracket algebra [26]. The study of projective 
geometry by Clifford algebra was initiated by Hestenes and Ziegler [18], and Stolfi 
[34]. The representation of projective transformations by spinors was initiated by 
Doran et al. [11], where a homomorphism of the Lie algebra gl[n ) into so{n,n) 
was proposed, making it possible to construct projective transformations by 
elements of Pin{n,n). Following this line, Goldman and Mann [15] discovered 
for many 3-D projective transformations their bivector generators in Cl( 4,4). 
Considering that the dimension of so(4,4) is Cf = 28, while the dimension of 
sl(4) is 15, the embedding space of sl( 4) seems too high [13]. 

A classical result [3] states that the group SL( 4), which acts upon the 4- 
space of homogeneous coordinates of points, is in fact isomorphic to the group 
Spino(3, 3), the connect component of Spin( 3,3) containing the identity, and the 
isomorphism is realized via the Pliicker coordinates of lines and the adjoint action 
of Spin( 3,3) upon R 3,3 . This canonical isomorphism indicates the possibility of 
using the wrench model, the model of spatial lines, to study projective geometry 
with Cl{ 3, 3). 

In AGACSE 2009, Li and Zhang [27] proposed a new model of 3-D projective 
geometry by taking the null vectors of R 3,3 as algebraic generators, and defining 
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points and planes as the two connected components of the set of null 3-spaces of 
R 3,3 respectively. Whenever an element of Spino(3,3) acts upon R 3,3 , it induces 
a projective transformation via the outermorphism of the action upon the null 
3-vectors representing 3-D points and planes. This approach was later followed 
by Klawitter [23], who proposed an explicit expression of the spinor inducing a 
projective transformation in 4 x 4 matrix form, and recently by Dorst [13], who 
constructed bivector generators for many 3-D projective transformations. 

When viewed from the homogeneous coordinates model R 4 of 3-D projective 
geometry, the R 3,3 model seems to have too many defects. The map from SL(4) 
to SO(3 ,3) is not surjective, nor injective. The projective transformations of 
negative determinant cannot be lifted to 0(3,3), and conversely, the elements 
of 0(3, 3) with negative determinant do not correspond to any projective trans¬ 
formation, but represent projective polarities where points are all mapped to 
planes. In the R 3,3 model, while lines are represented by vectors, the 3-D points 
and planes are represented by null 3-vectors, whose embedding vector space has 
dimension Of = 20. To make things worse, the mapping from R 4 to the null 3- 
vectors is quadratic, and defining the subset of null 3-vectors in the 20-D vector 
space they span is difficult. 

From the mathematical viewpoint, establishing the space R 4 of homogeneous 
coordinates from the 6-space R 3,3 spanned by lines requires rigorous mathemat¬ 
ical argument. It is the converse procedure of Pliicker’s construction of line co¬ 
ordinates from point coordinates. The well-definedness of the points and planes, 
and the covariance of the construction under suitable transformations of R 3,3 
need to be established. The benefits of using the null 3-vectors instead of the 
linear space R 4 to represent points need to be discovered. The groups S'Oo(3, 3) 
and Spino(3,3) are too small to cover the whole group of all 3-D projective 
transformations and polarities, and finding suitable covering groups to provide 
spin representations for all 3-D projective transformations and polarities is in¬ 
dispensable. 

So compared with other models of Geometric Algebra for classical geome¬ 
tries [17], [19], [24], [12], [26], [20], [14], [21], the line geometric model of 3-D 
projective geometry is much less developed. When every problem raised above 
is solved, then for the group of 3-D Euclidean transformations, a highly mature 
subject of study in Geometric Algebra, a one-to-one correspondence among the 
representations in the line geometric model and in other models need to be set 
up. 

As mentioned before, the screw algebra is valid only for Euclidean motions, 
and the corresponding group SE( 3) is only a subgroup of SL( 4). When SE(3) 
is replaced by another 6-D Lie subgroup of SL(4), then the Lie algebra se(3) 
of Euclidean motions is replaced by another 6-D Lie subalgebra of s/(4). Cor¬ 
respondingly, we can introduce new screw forms for the 6-D Lie subalgebra, 
together with the new “virtual work” of a wrench, which is still a vector of R 3,3 , 
along an infinitesimal “projective motion” represented by a screw form of the Lie 
subalgebra. The 6-D Lie subalgebra has its own Lie bracket, so the correspond¬ 
ing screw forms should have a different cross product. The new “virtual work” 
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should be related to the new cross product, or even be completely determined 
by it. 

We can go one step further by decomposing the 15-D algebra sZ(4) into the 
direct sum of five 3-D vector spaces, so that instead of using only pairs of 3- 
D vectors as in classical screw theory for the screw forms of se(3), we can use 
5-tuples of 3-D vectors to represent screw forms of sZ(4), and develop a “super¬ 
screw theory”, equipped with “super-cross product” and “super-virtual work”. 

For the purpose of developing a mathematically rigorous model out of the 
peculiar, unfamiliar and seemingly ineffective line geometric model, for more 
effectively describing and manipulating 3-D projective transformations with Ge¬ 
ometric Algebra, and with the ambition to further extend screw theory to pro¬ 
jective geometry, we picked up the research subject again in 2014, and after 
one-year’s hard work, we are confident to announce that from the algebraic 
viewpoint, this model is sufficiently mature now. The main contributions are 
summarized as follows: 

1. Rigorous establishment of the R 3,3 model for 3-D projective geometry. 

While in the classical model of projective line geometry only SL( 4) has spin 

representation, and all the spinors are in Spino(3,3), a rather unsatisfying lim¬ 
itation, the new model completely overcomes the limitation by providing pin 
group representations for all 3-D projective transformations and polarities, thus 
enlarging the transformation group four times. 

The group of linear regularities of R 3,3 is defined by 

RL( 3,3) := {B e GL( 3,3) | B T JB = ±J}, (1.3) 

where J is the matrix form of the metric of R 3,3 . Only when we computed 
the group acting upon the null 3-vectors induced by RL( 3,3) did we find the 
complete version of the line geometric model. The group RL{ 3,3) double covers 
the whole group of 3-D projective transformations and polarities in this manner, 
and the group Pin sp ( 3,3) quadruple-covers the latter, hence it can be used to 
construct versors for all kinds of 3-D projective transformations and polarities. 

The well-definedness of points and planes in the R 3,3 model, and the covari¬ 
ance of the representations are established. Some nice properties of reflections in 
R 3,3 are found, together with the classification of 3-D projective transformations 
induced by two reflections in R 3,3 . 

2. Construction of spinors in factored form inducing 3-D reflections and rigid- 
body motions, and discovery of the relation between the cross product of the 
screw forms of se(3) and the virtual work. 

For 3-D reflections and rigid-body motions, the spinors inducing them in 
Pin sp ( 3,3) in factored form are discovered. Since the bivector Lie algebra of 
Cl{ 3,3) is isomorphic to sZ(4), any element of se(3) has a bivector form, and 
the cross product of the bivectors equals the cross product of their screw forms 
as vectors of R 3,3 . On the other hand, a wrench is only a vector of R 3,3 . To 
make pairing with a bivector, a vector needs to be first upgraded to a bivector of 
T 2 (R 3 ' 3 ) by making inner product (tensor contraction) with a trivector. We show 
that this trivector is exactly the one complementary to the trivector defining the 
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cross product of the screw forms of se(3), and the latter trivector is exactly 
the lift of the quadratic form of R 3,0,1 . This correspondence shows the intrinsic 
the connection between the virtual work and the cross product of se(3). The 
connection between the wrench interpretation and the twist interpretation of 
the line geometric model is now clarified. 

3. Extension of the cross product and virtual work of the screw forms of se(3) 
to other 6 -D Lie subalgebras of sl(4). 

For many 6 -D Lie subalgebras of si (4), we have developed the corresponding 
screw forms together with the cross product and virtual work that are completely 
determined by the Lie bracket of the subalgebra. In particular, for so(K) where 
K is a quadratic form of R 4 with rank > 3, we have established the corresponding 
screw forms, cross products and virtual works, and discovered a striking fact: 
the trivectors for constructing new cross products and new virtual works are 
exactly the lifts of the quadratic form K by the Pliicker correspondence and the 
dual Pliicker correspondence to the trivector space. This result demonstrates 
that there is no intrinsic connection between the se(3)-interpretation and the 
wrench interpretation of line geometry. 

This paper is organized as follows. Section 2 is on the Pliicker model of 3-D 
projective geometry, the invariant group PR( 3) and its quadruple-covering group 
Pin sp ( 3, 3). Section 3 is on the covariance of the Pliicker transform and the dual 
Pliicker transform, and the explicit expression of the induced group element 
of PR( 3) by an element of RL( 3,3). Section 4 is on properties of projective 
transformations induced by one or two reflections of R 3,3 , and the construction of 
spinors in factored form inducing rigid-body motions. Section 5 is on the bivector 
representation of se(3), the screw forms and the cross products of screw forms. 
For reflections in space, the generating spinors are also constructed. Section 6 is 
on extension of classical screw theory to 6 -D Lie subalgebras of si (4) other than 
se(3). A framework of super-screw theory is also presented. 

2 The Pliicker model of 3-D projective geometry 

Definition 1 . A real n-space is said to have signature S, p ' q ’ r , if it has a basis 
with respect to which the metric of the n-space is diagonal, where the diagonal 
elements are composed of 1 of multiplicity p , and — 1 of multiplicity q , and 0 of 
multiplicity r = n — p — q. R p,? stands for R p,g ’°, called a non-degenerate inner- 
product space; R p stands for R p ’°’°, called a Euclidean inner-product space; R 0,0,r 
is called a null inner-product space, and R p>1 is called a Minkowski inner-product 
space. 

Definition 2. A Witt decomposition of R 3,3 refers to a decomposition of M 3,3 
into the direct sum of two null 3-spaces, say I 3 , J 3 . A Witt basis of the decom¬ 
position is composed of a basis of I 3 , say Ei,E 2 ,E 3 , and the corresponding 
Witt-dual basis of J 3 , denoted by E^, El,, Eg, such that E, • E' = Sij. 

Lemma 1. [8] R 3,3 has infinitely many Witt decompositions. Let R 3,3 = I 3 ©J 3 
be a fixed Witt decomposition o/M 3,3 . Then for any basis Ei, E 2 , E 3 of the 3-space 
I 3 , the corresponding Witt-dual basis E^E^Eg of the 3-space J 3 is unique. 
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Lemma 2 . Fix a Witt decomposition R 3,3 = I 3 ® J 3 . For any null 3-space S 3 
of R 3,3 , if the dimension n o/Sgfllg is even (or odd), then the dimension m of 
S 3 fl J 3 is odd (or even). 

Proof. When n = 0, then R 3,3 = I 3 ® S 3 . Let Ei, E 2 , E 3 be a basis of I 3 , and 
let the corresponding Witt-pairing bases in J 3 and S 3 be respectively E(, E 2 , E 3 
and si,S 2 ,S 3 . Then s» = E' + JT SijEj where Sij = —Sji. Let 


( 0 .Sl 2 Sl3\ 
—S12 0 S23 I 

—Sl3 ^23 0 ) 


If S = 0 then S 3 = J 3 and m = 3. If S 7 ^ 0, then its rank is 2, so its kernel has 
dimension 1. Let X = X 1 S 1 + Xi&i + X 3 S 3 G S 3 fl J 3 , then 0 = JT Xi(s t — E') = 
j so (xi,X 2 , x^) T is in the kernel of S T . This proves m = 1. 

When n = 2, let I3 be spanned by Ei,E 2 ,E 3 , and let S3 be spanned by 
Ei,E 2 ,s. Then s = S 3 E 3 + JTs'E(. Since S3 is a null 3-space, by s • Ei = 
s • E 2 = s 2 = 0 , we get s' x = S 2 = S3S3 = 0 . If S3 = 0 then S3 = I 3 , violating the 
assumption that n = 2. So S 3 7 ^ 0 and S 3 = 0 . We get s = S 3 E 3 and the 1-space 
spanned by E 3 is S3 fl J3. This proves m = 1. 

When n = 1 , let S 3 be spanned by Ei,s,t, such that s = S2E2 + S3E3 + 
s 2 E 2 + S 3 E 3 , and t = ^ 2 E 2 + f 3 Eg + t' 2 E 2 + t' 3 E' 3 . From s 2 = t 2 = s • t = 0, we 
get 

S2-S2 T S3S3 = t 2 t 2 d - ^ 3^3 = ^2^2 T -§3^3 -(- t 2 S 2 d - ^ 3^3 = 0 - ( 2 - 1 ) 


Since the 1-space spanned by Ei is I 3 fl S 3 , 


2 J 7 ^ 0 . 1 11 the 2 -space spanned 
*2 * 3 

1 \T 


by vectors (s 2 , s 3 ) T and (s' 2 , S 3 ) , by (2.1), we have (s 2 , S 3 ) T = A(sg, —s' 2 ) T , and 
(t 2 , ts) T = p{t 3 , —t' 2 ) T , and A = p. If A = 0 then S 3 is spanned by e^e^eg, so 


to = 2. If A 7 ^ 0, then 


s 2 S3 
t'2 t 3 


7 ^ 0, so S 3 fl J 3 = {0}, and m = 0. 


When n = 3, the conclusion m = 0 is trivial. 


□ 


Corollary 1. For a Witt decomposition R 3 ’ 3 = l 3 ©J 3 with Witt basis Ei, E 2 , E 3 
and E(, E 2 , E 3 , if S 3 fl I 3 is the 1-space spanned by E 1; then S 3 is spanned by 
Ei, s, t, where 

s — AE 2 T pE 3 , t = AE 3 — pE.^, (2.2) 

for some p 7 ^ 0 and A. If S 3 fl I 3 is the 2-space spanned by E 1; E 2 , then S 3 is 
spanned by E 1; E 2 , Eg. 

Lemma 3. Let the intersection of null 3-spaces P 3 .Q 3 be a 1-space, then for 
any null 3-space S 3 0 /R 3 ’ 3 , the dimension n of S 3 fl P 3 and the dimension m 
of S 3 n Q 3 have the same parity. 

Proof. Let P 3 fl Q 3 be the 1-space spanned by vector Ex. Assume S 3 7 ^ P 3 
and S 3 7 ^ Q 3 . 
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When n = 0, for the Witt decomposition R 3,3 = P 3 ® S 3 , by Lemma 2, m = 0 
or 2 . 

When n = 1 and P 3 (~l Q 3 = P 3 D S 3 , let Ei,E 2 ,E 3 be a basis of P 3 , 
and let E^E^Eg be the Witt-pairing basis of P 3 for a Witt decomposition 
R 3 ’ 3 = P 3 ffiP 3 . By Corollary 1, Q 3 is spanned by Ei, A 9 E 2 + ^ g E 3 , A 9 E 3 — /z^E^, 
where y q 7 ^ 0, while S 3 is spanned by Ei,A s E 2 + /z s E 3 ,A s E 3 — /^ S E 2 , where 
Us 7 ^ 0. Since S 3 7 ^ Q 3 , \ q : p q 7 ^ X s : y , s , so the 1-space E 3 is the only 
intersection of S 3 and Q 3 . This proves m — 1. 

When n = 1 but P 3 (~l Q 3 7 ^ P 3 fl S 3 , let Ei,E 2 ,E 3 be a basis of P 3 such 
that P 3 fl S 3 is the 1-space spanned by E 2 . Let E^E^Eg be the Witt-pairing 
basis of P 3 for a Witt decomposition R 3 ’ 3 = P 3 ffi P 3 . By Corollary 1, Q 3 is 
spanned by Ei, A g E 2 + ^ 9 E 3 , A g E 3 — /r g E 2 , where n q 7 ^ 0, while S 3 is spanned 
by E 2 , A s Ei + p s E' 3 , A s E 3 — /i.sEj, where /j s 7 ^ 0. Then S 3 fl Q 3 is the 1-space 
spanned by A s ^ 9 Ei + A 9 /x s E 2 + ^i s /i g E 3 . Again to = 1. 

When n = 2 and P 3 fl Q 3 C P 3 fl S 3 , let Ei, E 2 , E 3 be a basis of P 3 such 
that P 3 fl S 3 is spanned by Ei, E 2 . Let E(, E 2 , E 3 be the Witt-pairing basis of 
P 3 for a Witt decomposition R 3,3 = P 3 ffi P 3 , then S 3 is spanned by Ei, E 2 , E 3 , 
and Q 3 is spanned by Ei, AE 2 + /tE 3 , AE 3 — /jE( where fi 7 ^ 0. Then Q 3 fl S 3 is 
spanned by Ei, AE 2 + /jE 3 , and to = 2. 

When n = 2 but P 3 flQ 3 is not in P 3 fiS 3 , let Ei, E 2 , E 3 be a basis of P 3 such 
that P 3 n S 3 is spanned by Ei, E 2 , and P 3 D Q 3 is spanned by E 3 . Let E 3 , E 3 , E 3 
be the Witt-pairing basis of P 3 for a Witt decomposition R 3,3 = P 3 ffi P 3 . Then 
S 3 is spanned by Ei, E 2 , E 3 , and Q 3 is spanned by E 3 , AE 3 + /jE 2 , AE 2 — /iEj 

where n 7 ^ 0. Obviously Q 3 0 S 3 = {0}, and to = 0. EH 

Proposition 1. The set of null 3-spaces can be decomposed into two subsets 
(connected components): in each subset the dimension of the intersection sub¬ 
space of any two different elements is 1, and between the two subsets, the dimen¬ 
sion of the intersection subspace of any two elements, one from each subset, is 
0 or 2. 

As a corollary, the following concepts of points and planes in the set of null 
3-spaces of R 3,3 are well defined, and any pair of non-incident point and plane 
form a Witt decomposition of R 3,3 . 

Definition 3. For fixed Witt decomposition R 3,3 = I 3 ffi J 3 , if we call I 3 a point 
(the origin), and call J 3 a plane (the plane at infinity), then for any null 3-space 
S 3 of R 3,3 , let n be the dimension of the vector space S 3 fl I 3 , if n is even, then 
S 3 is called a plane, and if n is odd, then S 3 is called a point. 

The above concepts of points and planes have the following background. The 
classical Pliicker map changes a pair of points of Euclidean affine space £ 3 in their 
homogeneous coordinates form X = (xo, Xi, X 2 , xf) T and Y = {yo,yi,y 2 ,yf) T 
with respect to the basis eo, e 3 , e 2 , e 3 of R 4 , to a vector XY := XA Y of R 3,3 
where the induced basis is e l:] = e,e :] for 0 < i < j <3, and the image of the 
two points is represented by its Pliicker coordinates with respect to the induced 
basis. 
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Henceforth we always denote the outer product of the Grassmann algebra 
A(R 4 ) generated by R 4 by the juxtaposition of participating elements, while 
denoting the the outer product of the Grassmann algebra A(R 3,3 ) by the wedge 
symbol. 

Let eo,er,e 2 ,e 3 be an orthonormal basis of R 4 . Let 

Ei = e 0 i, E 2 = e 0 2 , E 3 = e 0 3 ; 

Ei = e 2 3 , E 2 = e 3 i, Eg = e- 12 . 

Then 

I 3 := Ei A E 2 A E 3 , J 3 := A Ej A E3 

are two null 3-spaces of R 3,3 forming a Witt decomposition. 

Let eo represent the origin of the Euclidean affine 3-space £ 3 , and let R 3 = 
(ei, e 2 , e 3 ) represent the plane at infinity. The basis eo, ei, e 2 , e 3 induces a Witt 
basis (2.3). The two 3-vectors I3, J3 are invariant under any special linear trans¬ 
formation of R 3 . 

For X, Y G R 4 , the vector XY G R 3,3 , if not zero, is a null vector. Conversely, 
any null vector of R 3,3 is the image of either an affine line or a line at infinity of 
£ 3 under the Pliicker map. A 2-space S 2 of R 3,3 spanned by null vectors, when 
interpreted geometrically so that its null 1 -spaces are lines in space, has two 
kinds: ( 1 ) a pair of non-intersection lines in space, when the signature of S 2 is 
R 1,1 ; (2) a pair of incident point and plane, i.e., a pencil of lines incident at a 
fixed point and at the same time lying on a fixed plane, when S 2 is null. A null 

3- space S3 of R 3,3 when interpreted in line geometry, represents either a point 
or a plane, i.e., either all lines incident at a fixed point, or all lines lying on a 
fixed plane. 

Definition 3 is based on a fixed basis e 0 , ei, e 2 , e 3 of R 4 and the induced Witt 
decomposition. Now Proposition 1 tells us that no matter what the underlying 

4- space R 4 is and what the induced Witt decomposition of R 3,3 could be, as 
long as the fixed null 3-space Ei A E 2 A E 3 of R 3,3 is classified as a “point”, 
then in a new line geometry whose abstract “lines” are the null 1-spaces of R 3 , 
any “point” defined with respect to the basis eo,ei,e 2 ,e 3 of the original R 4 is 
always classified as a point in the new line geometry. 

Fix the underlying space R 4 of the homogeneous coordinates of Euclidean 
affine space £ 3 , and fix a basis eo,ei,e 2 ,e 3 of it. The projective transformation 
group of £ 3 is the union SL( 4) U SL~( 4), where SL( 4) is the linear transforma¬ 
tions of determinant 1 , while SL~(4) is the linear transformations of determinant 
— 1. In the dual space (R 4 )* of R 4 equipped with the corresponding dual basis 
eg, e 3 , , e 3 such that the pairing between e.; and e* is Sij , the correspond¬ 
ing linear transformation of A G GL(4) is A~ T . The pair (A, A~ T ) acts upon 
R 4 x (R 4 )*, and is still called a general linear transformation. 

Any non-singular linear mapping from R 4 to (R 4 )* is called a projective 
polarity. The set of all projective polarities is denoted by GP(4). Such a mapping 
is called a special polarity if its matrix form A has determinant 1. The matrix 
A~ T represents the corresponding linear mapping from (R 4 )* to R 4 . The pair 
(A, A~ T ) acts upon R 4 x (R 4 )*, and is still called a projective polarity. The set 


(2.3) 

(2.4) 
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of special polarities is denoted by 5P(4), and the set of projective polarities with 
determinant —1 is denoted by SP~( 4). 

Obviously, the set 

UR( 4) := SL( 4) U SL~ (4) U SP{ 4) U SP~ (4) (2.5) 


is a group, and SL{4), SL(A)USL~(A), SL(4)USP(4) are three subgroups. XJR(A) 
is called the group of unitary regularities. 

On the other hand, consider some general linear transformations of R 3,3 . 
Fix a Witt decomposition I 3 © J 3 and the corresponding Witt basis Ei, E 2 , E 3 , 
E(, E(,, E 3 of R 3,3 . Let J be the linear transformation in R 3,3 interchanging E,; 
and E' for i = 1,2,3, i.e., its matrix form is 


J ■= 


( 0 I 3x3 \ 

\ ^3x3 0 ) ■ 


( 2 . 6 ) 


Let T be the linear transformation in R 3,3 changing E,; to — E* while preserving 
E' for * = 1,2,3, i.e., its matrix form is 


r := 



(2.7) 


Notice that T interchanges positive vectors and negative vectors of R 3,3 . The 
following is obvious: 

J r = -rj=(_° l3 I 0 3 ). ( 2 . 8 ) 

The group of non-singular linear transformations B in R 3,3 satisfying 

B T J r B = ±J (2.9) 

is called the group of linear regularities in R 3 ’ 3 , denoted by RL( 3,3). When 
the sign is positive, the corresponding subset forms the group of orthogonal 
transformations 0(3,3); when the sign is negative, the corresponding subset is 
called the anti-orthogonal transformations, denoted by AO(3,3). 

The subgroup of special orthogonal transformations SO( 3, 3) has two con¬ 
nected components: the component containing the identity transformation 16x6 
is denoted by 500(3,3), while the component containing — Ig X 6 is denoted by 
50r(3,3). 

J is an orthogonal transformations of determinant — 1. The subset of or¬ 
thogonal transformations of determinant —1 is denoted by SO~ (3,3). It also 
has two connected components: the component containing matrix J is denoted 
by 50()"(3, 3), while the component containing — J is denoted by 50f(3, 3). 

T is an orthogonal transformations of determinant — 1. The set of anti- 
orthogonal transformations of determinant —1 is denoted by SAO~{ 3,3). It has 
two connected components: the component containing matrix T is denoted by 
SAOq(3,3), while the component containing T is denoted by 3). 
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JT is an anti-orthogonal transformations of determinant 1. The set of anti- 
orthogonal transformations of determinant 1 is denoted by SAO(3,3). It has 
two connected components: the component containing matrix JT is denoted by 
S'AO 0 (3,3), while the component containing — JT is denoted by SAOi(3,3). 
We have 

RL( 3,3) = SO 0 (3,3) U SOi (3,3) U SO ^ (3,3) U SOf (3,3) , . 

U5AO 0 (3,3) U 5AOi(3,3) U SAOq (3,3) U SAOf(3, 3). ( Uj 

The two groups UR( 4) and RL( 3,3) are related by the Pliicker transform 
and dual Pliicker transform defined as follows. Let eo,ei,e 2 ,e 3 be a fixed basis 
of R 4 . Let zl 3 (R 4 ) be the realization space of (R 4 )*, whose basis 

eg = e! Ae 2 Ae 3 , ej; = -e 0 A e 2 A e 3 , ej = -e 0 Ae 3 Ae 1 , e 3 = -e 0 Ae 1 Ae 2 

( 2 . 11 ) 

satisfy for all positive permutations ijk of 123, the following: 

e 0 V ej = e* V e* = 1, 

e o V e* = e jk = E', (2.12) 

e* V e* = e ok = E fe . 

Here “V” is the meet product in the Grassmann-Cayley algebra generated by 
R 4 . The pairing between e,; and e* is defined by e p Ve*= 5 pq . Let (2.3) be the 
induced Witt basis of R 3,3 . 


Definition 4. The Pliicker transform from GL( 4) to GL( 3,3) is defined by A G 
GL{ 4) i-a A 2 A G GL( 3, 3), where 

(A 2 A)e,j = (Ae») A (Ae^) G R 3 ’ 3 . (2.13) 

The dual Pliicker transform from the set of projective polarities GD(4) to GL{ 3,3) 
is defined for any projective polarity D as 

(V 2 D)e. y - = (D ei ) V (De^) G R 3 ’ 3 . (2.14) 


For example, for the affine transformation 


A : 



Gi 4 G 




G R 4 , 


where t G R 3 and L G GL( 3), the matrix form of A 2 A with respect to the basis 
Ei, E 2 , E 3 , E ' l5 E' 2 , Eg is 

(txL L°0- (2 - 15) 

The following result is direct. 

Lemma 4. For any A G GL( 4), 


det(A 2 A) = (det(A)) 3 . 


(2.16) 
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For any D £ GD( 4), 

det(V 2 D) = —(det(D)) 3 . (2.17) 

Furthermore, let Ai, A 2 € GL{ 4) and Di,D 2 € GD(4), and let “o” denote the 
composition of mappings, then 

A 2 (Ai o A 2 ) = (A 2 Ai) o (A 2 A 2 ), V 2 (Ai o D 2 ) = (A 2 Ai) o (V 2 D 2 ), , . 

V 2 (D! o A 2 ) = (V 2 Di) o (A 2 A 2 ), A 2 (Dx o D 2 ) = (V 2 Dx) o (V 2 D 2 ). 

The Pliicker transform maps SL(4) onto SOo{3,3), and the kernel is ±I 4 X 4 - 
It also maps SL~( 4) onto SAOg(3, 3), such that ±A £ SL~( 4) are mapped 
to the same image of SAOff (3, 3). For example, the matrix diag(—1,1,1,1) is 
mapped to T £ SAOq(3, 3), and the pre-images of T are ±diag(— 1 , 1 , 1 , 1 ). The 
branches 5'Oi(3, 3) and (3,3) have no pre-image in SL{ 4) U SL~ (4). 

Similarly, the dual Pliicker transform maps SP~( 4) onto SOq (3,3), such 
that any B £ SOq (3, 3) has two pre-images ±A £ SP~( 4). For example, the 
mapping D : e; 1-4 e* for % = 0,1,2, 3 is mapped to J , and the pre-images of J 
are ±D. The dual Pliicker transform also maps <SP(4) onto SAOo(3,3), and the 
pre-images of any B £ SAOo(3,3) are of the form ±A £ SP( 4). The branches 
SAOi(3,3) and SOf( 3,3) have no pre-image in SP(4) U SP _ (4). 

Proposition 2. The Pliicker transform is a double-covering homomorphism from 
SL(4) to SOq( 3,3). The Pliicker transform and the dual Pliicker transform pro¬ 
vide a double-covering homomorphism in the sense of (2.18), from UR( 4) to the 
following subgroup of RL( 3,3): 

RL 0 ( 3,3) := SOq(3, 3) U SOq (3, 3) U S'AO 0 (3,3) U SAOq(3, 3). (2.19) 

In the setting of Clifford algebra Cl( 3,3), the Clifford product is always de¬ 
noted by juxtaposition of participating elements. Any element of the Pin group 
Pin( 3,3) is generated by invertible vectors of unit magnitude of R 3,3 under the 
Clifford product. Any element of the subgroup Spin( 3,3) is the Clifford prod¬ 
uct of even number of unit vectors. Spin~( 3,3) is the subset of elements that 
are the Clifford product of odd number of unit vectors. Spin( 3,3) has two con¬ 
nected components, and the component contain the identity element is denoted 
by Spino(3,3), the other component is denoted by Spini(3,3). 

Proposition 3. Any element of Spino(3, 3) is the Clifford product of even num¬ 
ber of negative vectors and even number of positive vectors; any element of 
Spini(3,3) is the Clifford product of odd number of negative vectors and odd 
number of positive vectors. In particular, ±1 are in Spino(3,3), while ± 13,3 are 
in Spini(3,3), where 

13,3 := Ei A E 2 A E 3 A Ei A El, A Eg = Ei 2 3 i/ 2 / 3 /, (2.20) 

is a pseduscalar o/yl(R 3 ’ 3 ) satisfying I 2 3 = 1. 
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Proof. Denote the set of positive vectors by P( 3,3), and denote the set of neg¬ 
ative vectors by N( 3, 3). First we prove that each of P(3, 3), JV(3, 3) is connected. 
By symmetry, we only consider P( 3, 3). Let vi, V 2 be two positive vectors, then 
vi A V 2 is one of R 2,0,0 , R 1,0,1 , R 1,1,0 . Hence there exists a positive vector V 3 that 
is orthogonal to both vi,V2- On Euclidean plane vi A V3, vi and V 3 are con¬ 
nected by positive vectors; on Euclidean plane V2 A V3, V2 and V3 are connected 
by positive vectors. So vi,v 2 are connected. 

Denote the set of elements that are the Clifford product of either two positive 
vectors or two negative vectors by P 2 ( 3 , 3), and denote the set of elements that 
are the Clifford product of either one negative vector and one positive vector, or 
one positive vector and one negative vector, by A^(3,3). Next we prove that each 
of P 2 (3,3), iV 2 (3,3) is connected. Let E + ,E_ be a pair of unit positive vector 
and unit negative vector in R 3,3 . Let 


v + =E + -AE_, v_ = E_ + AE_|_, (2-21) 

where 0 < A < 1. They are respectively a positive vector and a negative vector. 
By the continuity of the Clifford multiplication, we only need prove that E_v_ 
equals the Clifford product of two positive vectors, and E_E + equals the Clifford 
product of a positive vector and a negative vector. Both are true because 

E^v_ = -{E^(E_E+)}{(E+E_)(E_ + AE + )} = (-E+)(E+ - AE_), 
E_E + = (—E+)E_. 

( 2 . 22 ) 

Obviously ±1 G P 2 (3,3). By (E + E_)(E + E_) = 1 and the continuity of the 
Clifford multiplication, we get that the Clifford product of four invertible vectors, 
where the number of negative vectors is even, must be in the same connected 
component with P 2 ( 3 , 3). By (E + E_)(E + E + ) = E + E_ and the continuity of 
the Clifford multiplication, we get that the Clifford product of four invertible 
vectors, where the number of negative vectors is odd, must be in the same 
connected component with Afc)3,3). 

Denote by Spino(3,3) the connected component of Spin( 3,3) containing 
P 2 ( 3 , 3), and denote by Spin\(3 , 3) the connected component containing A^3, 3). 
By induction on the number of invertible vector factors in the factorization of 
an element of Spin(3,3), we get that Spino(3,3) contains all elements that are 
the Clifford product of even number of negative vectors and even number of 
positive vectors, while Spini(3,3) contains all elements that are the Clifford 
product of odd number of negative vectors and odd number of positive vectors. 
Since Spin( 3, 3) = Spino(3, 3) U Spin\{3, 3) and Spin( 3, 3) is known to have two 
connected components, Spino(3,3) and Spin\(3,3) are not the same connected 

component. EH 


Lemma 5. J € S'O - (3, 3) is double covered by ±(Ei — E' 1 )(Ei — E' 1 )(Ei — E^) G 
Pin~( 3,3), and —J is double covered by ±(Ei + E , 1 )(Ei + E' 1 )(Ei + E' x ) G 
Pin~( 3, 3). 
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Proposition 4. The set Spin~ (3,3) has two connected components. The com¬ 
ponent double-covering J is denoted by Spinff{ 3,3), whose elements are each the 
Clifford product of odd number of negative vectors and even number of positive 
vectors. The component double-covering —J is denoted by Spirif( 3, 3), whose el¬ 
ements are each the Clifford product of even number of negative vectors and odd 
number of positive vectors. In particular, all positive vectors are in Spinf (3,3), 
while all negative vectors are in Spinff (3,3). 

The group 0(3,3) is double-covered by Pin( 3,3), and the covering homo¬ 
morphism is given as following: for any ±U £ Pin( 3, 3), 

Ad{ jX := eUXU” 1 , for X £ R 3 ’ 3 (2.23) 

is a transformation belonging to 0(3,3), where e = 1 if U £ Spin(3,3), and 
e = —1 if U £ Spin~(3 , 3). 

Denote 

Pino(3 , 3) = Spino{3 , 3) U Spinff(3 , 3), Pin\{3, 3) = Spin\{3 ,3) U Spinf( 3,3). 

(2.24) 

Pmo(3, 3) is a subgroup of Pin( 3, 3), and double covers SOq{3, 3) U SOff{ 3,3). 
In particular, SOq{3,3 ) is double covered by Spino{3, 3). Now that the Pliicker 
transform provides another double-covering homomorphism of S'Oo(3, 3) by SL{ 4), 
we get the classical result that the two groups SL(4) and Spino(3,3) are iso¬ 
morphic. 

Below we extend the above double-covering map to Sb40(3,3) U SAO~ (3,3). 
We have seen that T £ <Sh40 _ (3,3). For any B £ 0(3,3), obviously TBT £ 
0(3,3), so any element of AO(3, 3) must have a unique matrix form TB for 
some B £ 0(3,3). In fact, we have the following: 

Definition 5. Define the following isomorphism in group Pin( 3,3): For any 
U = Yi Y 2 • • • Y r £ Pin( 3,3) where Y t £ R 3 ’ 3 , 

U r := (TY x )(TY 2 ) • • • (7Y r ). (2.25) 

Lemma 6 . For any U £ Pin( 3, 3) and any X £ R 3,3 , 

T(Ad{jX) = Ad* v r(rX). (2.26) 

Proof. If (2.26) is true for U = Yi, by induction it is true for any other 
element of Pin( 3, 3). It holds for U = Yi by direct verification. EH 

Definition 6 . By defining a formal associative product between T and Pin( 3,3) 
satisfying the following commutativity, a new group is generated, denoted by 
Pin sp ( 3, 3): for any U £ Pin( 3, 3), 

ToU = U r oT. (2.27) 

Let TPin( 3,3) be the coset of Pin( 3, 3) with respect to T, then 


Pin sp { 3,3) = Pin{ 3,3) U TPin( 3,3). 


(2.28) 
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The adjoint action of TPin( 3,3) upon R 3 ’ 3 is defined as follows: for any Uofo 
V G TPzn(3,3), where U, V G Pm(3,3), for any X G R 3,3 , 

A<j oToV X := Ad{j(T(Ad*vX)). (2.29) 

Pin sp ( 3, 3) double covers 0(3, 3)UAO(3,3) by the adjoint action. For i = 0,1, 
let TPinii 3, 3) be the coset of Pim( 3, 3) with respect to T for i — 0,1. Then 

Pin s 0 p {3 ,3) = Pm 0 (3,3) U TPin 0 ( 3,3) (2.30) 

double covers RLq(3, 3). Since UR( 4) also double covers RLq(3, 3) by the Pliicker 
transform and dual Pliicker transform, UR( 4) and PmQ P (3,3) are isomorphic. 
Below we realize T in Cl{ 3,3). Let 


K 2 = Eh/ + E 22 ' + E 33 / := eoi A e23 + eo 2 A e3i + eo3 A ei 2 - (2-31) 
It is called the symplectic form of R 3,3 with respect to the Witt decomposition 

R 3 ’ 3 = I 3 ® J 3 - 

Lemma 7. K 2 is invariant under any general linear transformation C in the 
3-space I3 and the associated linear transformation C~ T in the 3-space J3. In 
other words, it is independent of the the choice of Witt basis of the fixed Witt 
decomposition. 

Proof. Let ei, e 2 , e 3 , ej, e 2 , e' :i be a Witt basis of I3 ® J 3 , and let A G GL(l 3 ) 
such that Ae, = a.;. Let a,; = (au, a2i, a-3i) T , and let a' = {a' u , a' 2i , a' 3i ) T , where 
atj is the minor of A by removing the i-th row and j- th column. Then A~ T e( = 
a'/ det(A), and afa' = Sij det(A). We have 

ai A a' x ® ai A a 2 ® ai A a( 

= (onan' + 021021' + « 3 i^ 3 i') e i A e( + (anai2' + 021022' + « 3 i« 32 ') e i A e 2 

+(011013/ + 021023' + 03 i« 33 ') e i A eg + (012011' + «22«2i' + 032 « 3 i ') e 2 A e'i 

+ (012012' + 022 ^ 22 ' + 032032')e 2 A e 2 + (0l20l3' + « 22 « 23 ' + « 32 « 33') e 2 A eg 

+(ai3«ii' + a 23 a 2 i' + 033031')e3 A e( + (013012' + 023022' + 033032' )e3 A e 2 

+ (013013' + 023023' + 033033')eg A eg 

= det(A)(ei A e'g + e 2 A e 2 + e 3 A eg). [—j 

Some simple facts about K 2 : 

- K| = 3 - 2K 2 I 3 ,3- 

- K 2 is invertible: K ^ 1 = (K 2 + 2 I 3 i 3 )/ 3 . 

- For any X G R 3 - 3 , K 2 XK 2 = (1 -’2I 3 , 3 )X. 

- For any X, Y G R 3 ’ 3 , if X A K 2 = Y A K 2 , then X = Y. 

Proposition 5. For any X G R 3,3 , 


rx = x- K 2 . 


(2.32) 
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We have seen that UR{ 4) is isomorphic to only half of the group Pin sp ( 3,3). 
To study 3-D projective transformations and polarities, the homogeneous model 
R 4 of 3-D projective geometry does not allows for the whole group Pin sp ( 3, 3) to 
be used. To overcome this drawback, we need to proceed to use the null 3-space 
representation of projective points and planes, instead of returning to R 4 . This 
ideal leads to the following new model of 3-D projective geometry. 

Definition 7. The Pliicker model of 3-D projective geometry refers to the study 
of 3-D projective transformations and polarities by using the adjoint action of 
the whole group Pin sp ( 3, 3) upon the null 3-spaces of R 3,3 as the representation 
space of 3-D points and planes. For any element C 3 = Xi A X 2 A X 3 £ /1 3 (R 3,3 ) 
where X,; £ R 3,3 , for any U £ Pin sp (3,3), the adjoint action is the following: 

(A 3 Ad* v )C 3 := (Ad{ jXi) A (Ad{ jX 2 ) A (Ad{jX 3 ). (2.33) 

Definition 8. The 3-D projective regularity group PR(3) is defined as the quo¬ 
tient of UR( 4) modulo the equivalence relation in which two elements are 
equivalent if and only if they differ by scale. Usually we write 

PR{3) :=UR{A)/ * . (2.34) 

In the next section, it will be shown that via the Pliicker model, group 
Pin sp ( 3,3) quadruple-covers group PR( 3), and the kernel of the covering ho¬ 
momorphism is {±^± 13 ^}. 


3 Some properties of the Pliicker transform 

In this section we investigate two topics: null 3-vector representations of 3-D 
points and planes and their covariance under the group action of GL( 4) upon 
R 4 , and the adjoint action of Pin sp { 3,3) upon the null 3-vector representations. 

Notations of algebra. An orthonormal basis eo,ei,e 2 ,e 3 is fixed in R 4 , 
which induces a fixed Witt basis E t = eo* and E' = for positive permutation 
ijk of 123, where * = 1,2, 3. Denote 

Eq...^ =E 4l A---AE ir , (3.1) 

for ij £ {1, 2, 3,1', 2', 3'}. Denote ei 2 3 = eie 2 e 3 , and for any x £ R 3 , denote 

x 4- := x • ei 23 . (3.2) 

For example, e 4 - = ejk for positive permutation ijk of 123. 

For two 3-vectors P 3 , Q 3 £ /1(R 4 ), denote by P 3 Q 3 the meet product P 3 VQ 3 . 
This notation does not conflict with the same usage of the juxtaposition for 
representing the outer product of two vectors of R 4 , as the outer product of two 
3-vectors in yl(R 4 ) is always zero, so is the meet product of two vectors. 

Notations of geometry. Let x £ R 3 , then vector (xo,x) = xo e o + x £ R 4 
represents an affine point if and only if Xq ^ 0. A plane is determined by the 
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equation xq (—d) + n • x = 0 for point (xo, x) on it, where n G R 3 and d £ R, and 
at least one of them is nonzero, when dn ^ 0 , the equation represents the affine 
plane normal to vector n and with signed distance d/|n| from the origin along 
direction n. When d = 0, the plane passes through the origin; when n = 0, the 
plane is the plane at infinity. So the pair (n, —d) £ R 3 x represents the plane, 
and the representation is unique up to scale. 

Proposition 6 . Point (xo,x) £ R 3 has the following null 3-vector form in 

A 3 (R 3 ’ 3 ): 

(e 0 x) A ((e 0 x) • J 3 ) - x 0 (e 0 x) A K 2 + XqI 3 , ( 3 . 3 ) 

and plane (n, — d) has 3-vector form 

n 3- A (n 3 - • I 3 ) — dn 3 ^ A K 2 + d 2 J 3 , ( 3 . 4 ) 

Proof. When a;o e o + x is an affine point, where x = X\e\ + a; 2 e 2 + a; 3 e 3 , then 
Xq yf 0 , and the following three lines pass through the point, and form a basis of 
the 3-space of lines through the point: (xoeo + x)ei, (zoeo + x)e 2 , (xoeo + x)e 3 . 
We have 

((a; 0 e 0 + x)ei) A ((a; 0 e 0 + x)e 2 ) A ((a; 0 e 0 + x)e 3 ) 

= (zoeoi - x 2 e i2 + a; 3 e 3 i) A (x 0 e 02 + aqei2 - x 3 e 23 ) A (x 0 e 03 - Xie 3 i + x 2 e 23 ) 

= Z0I3 - ^ol^ieoi A (e 02 A e 3i + e 03 A e i2 ) + a: 2 e 02 A (e 0 i A e 23 + e 03 A e i2 ) 

+£3603 A (e 0 i A e 23 + e 02 A e 31 )} 

+®o( ;c i e oi + £2 e 02 + £ 3 e 03 ) A {(—aqe 12 A e 31 + x 2 e 12 A e 23 — x 3 e 3 i A e 23 } 

= a;gl 3 - a;§(e 0 x) A K 2 + x 0 (e 0 x) A {(e 0 x) • (e 23 A e 3i A e i2 )}. 

Removing the common factor Xo, we get (3.3). When Xo = 0, then e 0 x is the 
line through x and the origin, and (e 0 x) • J 3 are spanned by two lines at infinity 
that meet line eox, i.e., through the point at infinity x. So (eox) A ((eox) • J 3 ) 
is the 3-vector representing the point at infinity x. 

For an affine plane (n, —d) where n 2 = 1, let a, b, n be an orthonormal frame 
of R 3 . Then rfo = ab, and 

I 3 = (e 0 n) A (e 0 a) A (e 0 b), 

J 3 = (na) A (nb) A (ab), 

K 2 = (e 0 n) A (ab) - (e 0 a) A (nb) + (e 0 b) A (na). 

The following three lines span the 3-space of lines on the plane: (e<j + dn)a, (eo + 
dn)b, ab. We have 

((e 0 + dn)a) A ((e 0 + dn)b) A (ab) 

= (e 0 a) A (e 0 b) A (ab) + d 2 (na) A (nb) A (ab) 

+d(ab) A {(eoa) A (nb) — (e 0 b) A (na)} 

= (ab) A ((ab) • I 3 ) - d(ab) A K 2 + d 2 J 3 . 

Thus we get (3.4) for affine plane. When n = 0, then J 3 is the plane at infinity. 

□ 
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Definition 9. The map 

T : a; 0 e 0 + x (e 0 x) A ((e 0 x) • J 3 ) - a; 0 (e 0 x) A K 2 + *§I 3 (3.5) 

is called the Pliicker representation of the point. The plane (n, d) can be repre¬ 
sented by vector deo + n; for the plane, the map 

T' : de 0 + n i-*- n^ A (n^ • I 3 ) + dn x A K 2 + d 2 J 3 (3.6) 

is called the Pliicker representation of the plane. 

Remark. The Pliicker representations of points and points at infinity eo,ei, 
e 2 ,e 3 are respectively Ei 23 ,Ei 2 / 3 /,E 23 /i/,E 3 i/ 2 /. Furthermore, 

F{x 0 e 0 + x) = ^(croeo) + J'(x) - ((cr 0 e 0 )x) A K 2 , 

T'{de o + n) = F'(de 0 ) + ^'(n) + (dir 1 ) A K 2 . 


Definition 9 relies upon the specific Witt decomposition I 3 ©J 3 . The following 
property gives an alternative definition of T that is independent of the Witt 
decomposition. 

Proposition 7. Let X, a 3 , a 2 , a 3 be a basis o/M 3 . Then 


t-(- y-\ _ (Xa 3 ) A (Xa 2 ) A (Xa 3 ) 
1 [Xa ia 2 a 3 ] 


(3.8) 


In particular, .F(X) is independent of the choice o/ai,a 2 ,a 3 as long as the latter 
forms a basis of the 4~space together with X. Furthermore, 


T' = (A 3 J) o T. 


(3.9) 


Proof. (3.9) is obvious by (3.5) and (3.6). We only need to prove (3.8). When 
X is an affine point, by definition, 

(Xe 3 ) A (Xe 2 ) A (Xe 3 ) 

1 J [Xe ie 2 e 3 ] 

For * = 1,2,3, let a.^ = d;X + Ci,ei + c 2 ie 2 + c 3 ie 3 . Then det(cjj)j i j = i .. 3 ^ 0, as 

'1 d T ' 


(X ai a 2 a 3 ) = (X e 3 e 2 e 3 ) 


0 C 


where d T = (di,d 2 ,d 3 ), and C = {cjf)ij- 1 .. 3 . By direct verification, 

(XaQ A (Xa 2 ) A (Xa 3 ) _ det(C)(X ei ) A (Xe 2 ) A (Xe 3 ) 

[Xaia 2 a 3 ] clet(C)[Xeie 2 e 3 ] 

When X is a point at infinity, let X = Ax, where x is a unit vector, and 
let x, a, b be an orthonormal basis of the 3-space spanned by ei,e 2 ,e 3 , then 
[Xe 0 ab] = —A. By definition, 


T'(X) = A 2 (e 0 x) A ((eox) • (ab A bx A xa)) 


(Xe 0 ) A (Xb) A (Xa) 

[Xeoba] HH 


The following proposition answers the following question: Given J^X) for 
X = soeo + x, how to recover X? 
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Proposition 8 . For any X = a;o e o + x G R 4 , and any 77 = de q + n* G (R 4 )*, 

J r (X)-K 2 = - 2 a; 0 (e 0 x), 

IF" (II) • K 2 = -2d (e* Q n*). 1 


Proof. Let x = xiei + X 2&2 + ^ 3 ^ 3 , then 

3 3 

F(X) ■ K 2 = ^ Xi(Ei A (e 0 x)) • J 3 - 3x 0 (e 0 x) + x 0 = -2cc 0 (e 0 x). 

The second result can be proved similarly. EH 

The following is on the covariance of .F(X): if X G R 4 undergoes a general 
linear transformation, how does FIX) change in yl 3 (R 3,3 ) accordingly? 

Proposition 9. Let A be a non-singular linear transformation in the 4~space 
spanned by eo, ei, e 2 , e 3 . Then for any vector X = Xoeo + x of the f-space, 

det(A).F(AX) = A 3 (A 2 A).F(X), (3.11) 

and 

det(A)F'(A~ T X) = A 3 {J o (A 2 A" t ) o J)F'(X). (3.12) 

Proof. Without loss of generality, let X represent an affine point in space, 
then 


A 3 (A 2 A).F(X) 


(A 2 A)(X ei ) A (A 2 A)(Xe 2 ) A (A 2 A)(Xe 3 ) 


det(A) 


[Xeie 2 e 3 ] 

(AX)(AeO A (AX)(Ae 2 ) A (AX)(Ae 3 ) 


[(AX)(A ei )(Ae 2 )(Ae 3 )] 


det(A)J r (AX). 


(3.12) is obvious from (3.11), where A T is the induced linear transformation 
in the vector space of planes. EH 

Alternatively, a plane (n, d) can be represented by a 3-vector in A(R 4 ). Let 
n = aei + &e 2 + ce 3 , then 77 = deg + ae^ + foe 2 + ceg represents the plane. A 
point X = a"o e o + a’iei + a: 2 e 2 + a: 3 e 3 is on the plane if and only if X V 77 = 0, 
i.e., xod + x\a + a; 2 & + ai 3 c = 0 . 

Proposition 10. Let A G GL(A) act upon R 4 . In yl 3 (R 4 ) and yl 2 (R 3 ’ 3 ), the two 
matrices A 3 A = det(A)A~ T and B = V 2 (A 3 A) satisfy 

det(B) = (det(det(A)A ~ T )) 3 =det(A) 9 . (3.13) 

Let n* = aej + be?> + ceg. When 77 = deg + n* is an affine plane not through 
the origin, then dn* ^ 0 , and the plane meets the following three planes passing 
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through the origin: e|,ej,e|, so lines lie* for i = 1,2,3 form a basis of the 
3-space of lines on plane II. By (2.12), 

{net) A (ne* 2 ) A (77e§) 

= (de 23 — &eo 3 + ceo 2 ) A (de 3 i + aeo 3 — ceoi) A (de i 2 — ae q 2 + be or) 

= d 3 J 3 + c? 2 {ae23 A (eo 2 A e 3 i + eo 3 A ei 2 ) + be 31 A (eoi A e23 + eo 3 A ei 2 ) 

+cei 2 A (eoi A e23 + eo 2 A e 3 i)} 

+d(ae23 + foe 3 i + cei 2 ) A {(aeo 2 A eo 3 + 6 eo 3 A eoi + ceoi A 602 } 

= d 3 J 3 + d 2 (e* 0 n) A K 2 + d(e* o n) A {(eJTT) • I 3 }. 

Removing the common factor d, we get the following alternative form of (3.4): 

[e*n) A ((e*77) • I 3 ) + d(e*II) A K 2 + d 2 J 3 . (3.14) 

When 77 passes through the origin, then d = 0, and (3.14) still represents the 
plane. When II is the plane at infinity, then II = dej, so eJTI = 0, and (3.14) 
once again represents the plane. 

Definition 10. When a plane (n, d) is represented by a 3-vector 77 = de 3 + n*, 
the map 

T" :de*+n* ^ (e*II) A ((e^TI) • I 3 ) + d{e*II) A K 2 + d 2 J 3 (3.15) 
is called the dual Plucker representation of the plane. 

Remark. 


T"(de* 0 + n*) = T"{de*) + JT"(n*) + ((de5)n*) A K 2 . (3.16) 

Lemma 8. For the plane (n, d) represented by 3-vector 77 = dej) + n*, 

d = ne{e* 2 e* 3 = 77 V e* V e^ V e*. (3.17) 

Lemma 9. Let A £ GL(4 ) act upon M 4 . Then in R 3,3 , 

A 2 A = det(A) (V 2 A- t ). (3.18) 

Proof. Let A = (a 0 ai a 2 a 3 ). On one hand, A 2 A(e 0 i) = A(eo) A A(e*) = 
ao A a i, and A 2 A(e^) = a, A a.j. On the other hand, for positive permutation 
ijk of 123, 

V 2 A~ T (e 0i ) = V 2 (det(A) _1 A 3 A)(e* V e£) 

= det(A) _2 (a 0 A a*, A a,) V (a 0 A a; A a j) 

= det(A) _1 a 0 A a^; 

V 2 A- T (ey) = V 2 (det(A) _1 A 3 A)(e* V e*) 

= — det(A) _2 (a; A a, A a*,) V (a 0 A a, A a,) 

= det(A) _1 a, A aj. Q 
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Proposition 11. For Xel 4 and U G A 3 (R 4 ), 

.F(X) • F"{II) = -[XU] 2 . (3.19) 


Proof. Let X = Xoeo + x and U = deg + n*, then 

■F(X) • T"(n) = (J-(X) + x 2 J-(eo) - x 0 (e 0 x) A K 2 ) 

•(^"(n*) + d 2 T"{e* 0 ) + d(e* 0 n*) A K 2 ) 

= — (x • n) 2 — x^d 2 — 2 xq d (x • n) 

= - [XU] 2 . □ 


The proof of the following proposition is similar to those of Proposition 7 
and Proposition 9, and is omitted. 

Proposition 12. Let U, aj, a?i, a); be a basis o/(R 4 )*. Then 


F'\II) 


(Ua*) A (Ua*) A (Ua|) 
Ua]a*^ 


(3.20) 


In particular, F"(II) is independent of the choice o/a*,aJ,a 3 as /on <7 as f/ie 
latter forms a basis of the f-space together with U. 

Let A be a non-singular linear transformation in R 4 . It induces a linear 
transformation A~ T = def(A) -1 A 3 A in (R 4 )*. For any LI = de g + n* £ (R 4 )*, 

det(A^ T )J r ' / (A~ T U) = A 3 (V 2 A~ T ).F"(U). (3.21) 


Below we extend the previous results from GL{ 4) to GU(4). The following is 
a simple fact from linear algebra on the influence of coordinate transformations 
in R 4 upon the matrix form of A : R 4 —> (R 4 )*. 

Proposition 13. Let C = (co Ci c 2 C 3 ) £ GL( 4) be a coordinate transformation 
in R 4 , and let A be an invertible linear mapping from R 4 to (R 4 )*, whose matrix 
form with respect to the basis e.; ’s 0 /R 4 and e* ’s of (R 4 )* is A± X i. Then with 
respect to the basis ci’s of R 4 and its dual basis c*’s of (R 4 )*, where c* = 
(A 3 C)e*/det(C), the matrix form of A is C T A 4 X 4 C. 

Proof. The coordinate transformation C in R 4 induces the following coordi¬ 
nate transformation in yl 3 (R 4 ): 

3 

(A 3 C)e* = det(C)(C- T ) jie *. (3.22) 

i=o 

Let Ae* = J2j=o a ij e h and let Ci = ( c o*> c u, c 2 i, c 3 *) T . Then 
33 3 

Ac i — ^ ^ c^Ag/j — ^ ^ C’kiQ'jk^j — det(C) - 1 E Cki Qj k (C T )y(A 3 C)er 

k—0 j,k —0 j,k,l—0 

3 

= Y,(C T AC)ucl 

1=0 
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□ 


Notation. For any A 4 X 4 = (a>ij)ij,= 0 .. 3 ? denote its 2 x 2 minor composed of 
the u,v rows and p,q columns by AJJ^, i.e 


A uv 

pq * 


dup O'uq 
Clyp CLyq 


(3.23) 


Let A £ GD{ 4 ), Le., a non-singular linear map from R 4 to (R 4 )*. Let the 
matrix form of A with respect to the basis eo, ei, e 2 , e 3 of R 4 and the dual basis 
e o, e|, ej, eg of (R 4 )* be A 4x4 = (a 0 a x a 2 a 3 ), where a , = (ao„«i„o 2i ,«3i) T . 
Then Ac, = a* := J2j=o (l .n e j- and 

(V 2 A)(e,j) = a* a* = £ Af/e rs , (3.24) 

0<p<q<3 


where pqrs is a positive permutation of 0123. 


Lemma 10 . For the 4x4 identity matrix I 4x 4 representing an element of 
GD( 4), V 2 I 4x4 = J. 

Proposition 14. Let A £ GD( 4) have matrix form A 4x4 , and Let B £ GL( 3,3) 
be the matrix form o/V 2 A. When matrix A 4x4 is taken as a linear transforma¬ 
tion of R 4 , then A 2 A 4x4 has the matrix form J B £ GL( 3,3). 

Proof. Let A be the non-singular linear transformation in R 4 having the same 
matrix form A 4x4 = (ao a 4 a 2 ag) with A. Then the matrix form B of A 2 A is 
given by 

(A 2 A)( eii )= £ A%e pq . (3.25) 

0<p<q<3 

Obviously B = J'B. EH 

Proposition 15. Let A £ GD( 4 ). Then for any X £ R 4 , 

det(A)J r "(AX) = A 3 (V 2 A)J-(X). (3.26) 


Proof. We have 


A 3 (V 2 A)J-(X) 


(AX)(Aei) A (AX)(Ae 2 ) A (AX)(Ae 3 ) 
[Xeie 2 e 3 ] 


Since (Ae 0 )(Ae!)(Ae 2 )(Aeg) = aQa*a 2 a| = det(A)eQege5eg = det(A), we 
have det(A)[Xeie 2 e 3 ] = (AX)(Aei)(Ae 2 )(Ae 3 ). On the other hand, 


= (AXXAei) A (AX)(Ae 2 ) A (AX)(Ae 3 ) 

1 j (AX) (Ae 4 ) (Ae 2 ) (Ae 3 ) ‘ □ 


For any U £ Pin sp ( 3,3), A 3 Ad{j maps a null vector to a null vector. If 
U £ Spin{ 3,3) UTS'pm(3,3), then A 3 Ad\j maps the set of null 3-spaces repre¬ 
senting points to the same set; if U £ Spin~(3, 3) U TSpin~(3, 3), then A 3 Aci;jj 
interchanges the set of null 3-spaces representing points and the set of null 3- 
spaces representing planes. 
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Lemma 11. For any X £ R 4 and II £ (R 4 )*, 

.F(X) • I 3 ,3 = -F(x), r\n) ■ I 3 ,3 = -F"(n). ( 3 . 27 ) 

Furthermore, for any P £ R 3,3 , Ad^ 3 3 P = — P. 

Given a null 3-vector, (3.27) can be used to determine whether it represents 
a point or plane of £ 3 . 

Notation. For two algebraic elements a, b, if they differ by scale, we write 

a = b. (3.28) 

In ff 3 (R 3 ’ 3 ), the representation of a point or plane by a null 3-vector is unique 
up to scale. So Ad\ ( leaves each point and plane of £ 3 invariant, and induces the 
projective transformation represented by the the identity transformation I 4 X 4 of 
R 4 . Asa consequence, the adjoint action of Pin sp ( 3,3) upon the null 3-vectors of 
A 3 (R 3 ’ 3 ) induces a quadruple-covering homomorphism upon the group PP(3), 
and the kernel is composed of 4 elements: ±1, ±I 3j3 . 

Definition 11. For U £ Pin sp ( 3,3), The 4x4 matrix representation of Ad {j 
refers to a matrix form of the linear map induced by Ad\j via the null 3-vector 
representation of points and planes, from R 4 to R 4 or (R 4 )*. The matrix form is 
unique up to scale. 

As Pin sp ( 3,3) double covers RL( 3,3), any element B £ RL{ 3,3) also has 
4x4 matrix representation. By the covariance of T and T ", the 4x4 matrix 
representation after rescaling is mapped to one of ±B by the Pliicker transform 
or dual Pliicker transform. We study how to compute the 4x4 matrix form 
A of B. Without loss of generality, let B £ S'Oo(3, 3). Then B is the image of 
A £ SL( 4) under the Pliicker transform. 

Notations used only in this section. Denote 

P(0) := 123, P(l) := 12'3', P( 2) := 23'T, P(3) := 3l'2'. (3.29) 

Let ijk be a positive permutation of 123, denote 

Q(i) ■- j'k', Q{ij) := k'. (3.30) 

Proposition 16. Let A = (ai a 2 a 3 a 4 ) = (a pg ) p>g= o ..3 G SL(4), and let B = 
(^a/ 3 )a,/ 3 =i, 2 , 3 ,i', 2 , , 3 ' be the matrix of A 2 A with respect to the basis Ei,E 2 ,E 3 , 
Ei/, E 2 /, E 3 /, i.e., for 1 < i,j < 3, let ipq and juv be both positive permutations 
of 123, then 

h — A 0i h-, — A pg — A 0i \pi 13 311 

Denote by B the 3x3 minor formed by the u, v, w rows and p, q , r columns 
of B. Then for any 0 <i,j <3, 


det(A)o?- = B 


(3.32) 
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Proof. Without loss of generality, consider j 
commutes: 


ei 

1 


Ei 2 / 3 / 

1 

A| 


| a 3 b 

1 


1 


det(A) J 7 


ai 

- > 

C 3 


where 


1. The following diagram 


(3.33) 


C 3 = det(A){aoil3 - a 0 i(e 0 ai) A K 2 + (e 0 ai) A ((e 0 ai) • J 3 )} 


— det(A){aQ 1 Ei 2 3 + aoi&n(E 3 i 3 ' — E122') + aoi« 2 i(Ei 2 i' — E 233 ') 
+aoia 3 i(E 232 ' — E 3 h') + OiiEi 2 ' 3 ' + 02 1 E 23 'i' + a 3 iE 3 i' 2 ' 

+ana 2 i(Ei 3 'i' + E 22 ' 3 ') + aua 3 i(Eu< 2 ' + E 32 ' 3 ') + 021031 ^21/2' + E 33 /]/)}. 

(3.34) 

Since A 3 B(Ei 2 / 3 /) = J 2 ijk ®i2'3'E»jfc, by comparing the coefficients of Ei 23 , 
Ei 2 ' 3 ', E 23 'i', E 3 i/ 2 ' in this expression and in C 3 , we get 

= det(A)ag 1 , = det(A)of 1 , 

B^i; =det(A)ali, = det(A)o§ 1 . □ 


Remark. By (3.32), a matrix B G SO( 3,3) is in 500(3,3) if and only if 
Bpj*) > 0 for all 0 < i, j < 3. 

Proposition 17. Let ijk be a positive permutation of 123, and let l be one of 
0,1,2,3, then 

det{A)aoiau = -B^ = -B^', 

det(A) W = B’^if = B«j;W; 
det(A)a io azi = -bJJ, } = -B%$, 
det(A )a liaij = Bf g W= B J«.. 

Proof. Without loss of generality, consider l = 1. In the diagram (3.33), by 
comparing the coefficients of E i22 /, E 3 i 3 /, Ei 2 i-, E 233 /, E 3 n-, E 232 / in J] B*/ 2 , 3 ,E ijk 
and C 3 of (3.34), we get the first line of (3.35) for l = 1. By comparing the co¬ 
efficients of Ei 3 /i/,E 22 / 3 /,En/ 2 /,E 32 / 3 /,E 2 i/ 2 /,E 33 /i/, we get the second line of 
(3.35) for l = 1. When replacing A with A 1 , by the trivial fact 

A 2 (A t ) = (A 2 A) t , (3.36) 


we get the last two lines of (3.35). 


□ 


Proposition 18. Let ijk and pqr be both positive permutations of 123, then 
det(A )(ao p (ii q + a 0q a ip ) = —B^^ — B^^ = — B*q^ — B g ^ ? , 

(3.3 

AatfWfr, n I rt rt \ TD^QO") I TD^QU) p>Q(*)j 1 

aet(A){ai P ajq + a.i q aj P ) — B pg ^ + B g ^ g — ** p Q( q ) + **Q( P ) q - 
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Proof. Without loss of generality, consider i = 1 and j = 2. The following 
diagram commutes: 


ei + e 2 


Ei2'3' + E 2 3'i' + E 13 / 1 / + E 22 '3' 


A | |A 3 B 

i i 


ai + a 2 


det(A).F 


D 3 


(3.38) 


where for C 3 of (3.34), 


D 3 — C 3 + det(A){ 2 aoiao 2 Ei 2 3 + (aoiai 2 + ao2<iii)(E3i 3 / — Ei 22 /) 

+(a 0 ia 2 2 + a 02 a 2 i)(E 121 / — E 233 /) + (a 0 ia 32 + ao 2 a 3 i)(E 232 / — E 3 n/) 

+2anai 2 E 12 / 3 / + 2a 2 ia 22 E 2 3/i/ + 2a3ia3 2 E 31 / 2 / 

+(ana 2 2 + ai 2 « 2 i)(Ei 3 'i' + E 22 ' 3 ') + (ana 32 + 012031) (En/ 2 ' + £32/3/) 

+(021032 + 0,220,31) (E2V2 1 + E33/1/)}. 

(3.39) 

Since A 3 B(E 13 /;p + E 22 ' 3 ') = J2ijk( B i3'i' + B 22'3 ') E yfc, b y comparing the 
coefficients of E 3 i 3 >, Ei 22 ', Ei 2 i', E 233 ', E 232 ', E 311 /, E 13 / 1 /, E 22 ' 3 '> En/ 2 ', E 3 2 ' 3 ', 
E 2 i' 2 ', E 33 / 1 / in this expression and in D 3 — C 3 , we get for positive permutation 
pqr of 123 the following: 


det(A)(a 0 io p2 + o 02 a P i) = -B™*, - B™£, = -B& v - B^,, 
det(A)(o pl o ?2 + o 92 o p i) = B^;?; + B£# = B&* + B&*. 

□ 

The 4x4 matrix representation A 4 X 4 = (o pq )p >g = 0..3 of B S SOo(3,3) is 
given explicitly by (3.32), (3.35), and (3.37) as following: one entry of the first 
row of A must be nonzero, say ao qo 0 for some 0 < go < 3. From (3.32) we get 
Oq 9o ; from ao qo and (3.35) we get a,; 9o for all 1 < i < 3, and oo 9 for all q y^ go- 
From ao qo Oi q + ai qo cio q of (3.37) we get the rest components Oi q of A. The last 
step can be made more explicit as follows. 

By 


OQpOiq OQqOip — A p q — bi r > , 
OipCljq OiqOjp — A p q — bp / r ', 


we get 
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Proposition 19. 


ZaopCiiq — B p Q( g ) B Q^ q + biQ(pq) det(A) 

= -B& - B %« )q + b iQ{pq) det(A), 

2ao q a ip = —B — B g^ q — biQ( pq ) det(A) 

= - B pQM Bjg), - b lQ(pq) det(A), 

2 ai P aj q = BpQ*^) + B+ ^Q(y)Q( P g) det(A) 
= B pQ( 9 ) + B Q(p)g + b Q(ij)Q(pq ) det(A), 
2 iCLiqCljp = B pQW ^Q(p)q bQ(ij)Q( pq ) det(A) 

= Bg^j + - 6 Qfe)Q(p9) det(A). 


4 From reflections in M 3,3 to rigid-body motions in £ 3 

Any non-null vector X € R 3,3 generate a reflection of R 3,3 by its adjoint action. 
The induced action in the null 3-spaces interchanges points and planes, so it is 
a 3-D projective polarity from R 4 to (R 4 )*. 

Notations. For any x = ( x,y,z ) £ R 3 , denote 

E(x) := xEi\ + yE 2 + 2 E 3 , (ai\ 

E'(x) := xE[ + yE ' 2 + zE' 3 . 1 ' 

Definition 12. Any vector X £ R 3,3 must be of the form E(x) + E'(y) for 
x, y £ R 3 . The vector 

= (x, y) T £ I 3 x I 3 (4.2) 

is called the screw form of vector X. 

One advantage of the screw form is its capability of representing two different 
inner products simultaneously by the same dot symbol with no ambiguity. On 
one hand, M 3 is naturally equipped with the standard Euclidean inner product. 
On the other hand, the inner product of two vectors Xi = E(xi) + E'(yi) 
and X 2 = E(x 2 ) + E'(y 2 ) of R 3,3 is denoted by the same dot symbol. The two 
inner products are distinguished by the participating vectors, with the lower-case 
letters denoting vectors of R 3 , while the capitals denoting vectors of R 3,3 . 

We have 

Xl ' X2= (yl) ' (y 2 ) = xi -y»+yi- x 2- 

Vector Xi is invertible if and only if xi • yi ^ 0. In 4(R 4 ), E'(x) is the dual of 
E(x). In R 3 ’ 3 , 



E(x) • E'(y) = x • y. 


(4.4) 
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Notations. For two vectors x, y £ R 3 , x-y and xx y are the usual Euclidean 
inner product and vector cross product in vector algebra. For a 3 x 3 matrix 
M = (nil m 2 m 3 ) where each irq £ R 3 , denote 


x x M := (x x nil x x m 2 x x m 3 ). 
Alternatively, matrix x x M is defined by 

(x x M)y = x x (My), for y £ R 3 . 
For example, when x = (x , y, z ) T , 


0 


x x I 3 x3 = 


0 —x 


(4.5) 


(4.6) 


(4.7) 


Let U = (x, y) T £ R 3,3 be invertible, then for any V = (p, q) T £ 


V U 

Ad* u (V)=V-2 w - JJ U = 


/ xq + py \ 

/ P- x \ 

x-y 


xq + py 

q-y 

x-y 


By direct computation, we get 

Proposition 20. For invertible U = (x,y) T £ R 3,3 , 


(A 3 A^)J 3 = (/\ 3 Ad{j) J-(x) = -(x • y)J 3 . 

x-y 


(4.8) 


(4.9) 


In other words, the plane at infinity is mapped to the point at infinity x by the 
adjoint action of U. Furthermore, with respect to the basis e; of R 4 and e* of 
(R 4 )*, the 4x4 matrix form of Ad{j of unit determinant is 


±- 


1 


0 —y T 

y xx l 3x3 


\/l x -y| 

The matrix is skew-symmetric. Its inverse-transpose is 

1 


±- 


\/l x 'yl 


0 -x T 
X y X I 3x3 


(4.10) 


(4.11) 


Ad\j induces an affine transformation of £ 3 if and only if it preserves the 
plane at infinity. In yl(R 3,3 ), the plane at infinity is represented by J 3 = E]/ 2 ' 3 '. 
The following is a direct corollary of the above proposition. 

Corollary 2. Let U,; = (x.^, y,) T £ R 3 ’ 3 for i = 1,2 be invertible. Then ^4duiU 2 
induces an affine transformation of £ 3 if and only if Xj x X 2 = 0. If Xi = X 2 = x, 
then the affine transformation is volume-preserving if and only if x • y! = x • y 2 . 
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Definition 13. Any element M £ /1 2 (R 4 ) = R 3,3 induces a linear transforma¬ 
tion of R 4 as following: 


£(M) : XeR 4 £M'XeI 4 , (4.12) 

where the inner product is in C1(R 4 ). It is called the bilinear form of vector 
M £ R 3 ’ 3 . 


The matrix form of the transformation £(M) is obviously skew-symmetric. 
The following lemma is by direct verification: 

Lemma 12. For any M = (x, y) T G A 2 (R 4 ), 


C 



0 -x T \ 
x y x I 3x3 J ’ 


(4.13) 


When M £ R 3,3 is invertible, the above matrix is exactly the 4x4 matrix of 
Ad^i from the space of planes (R 4 )* to the space of points R 4 . 

The inner product in the space of 4 x 4 matrices is that of the embedded 
space R 16 of components. For any Mi, M 2 £ R 3,3 , 

£(M0 • C(JM 2 ) := tr ((£(M 1 )) T (£( l /M 2 ))) = 2(M 1 ■ M 2 ). (4.14) 

Corollary 3. Let M = (x, y) T £ R 3,3 . When x • y / 0, then ker(£(M)) = 0; 
when x • y = 0 but x, y are not both zero, then ker(£(M)) is the 2-space of R 4 
containing all the projective points of the line with Pliicker coordinates (x, y) T . 

An element of Spin(3 , 3) in factored form with respect to the Clifford product 
has the following nice property: 

Proposition 21. Let U = ViV 2 ---V r £ Pin( 3,3), where each V, £ R 3,3 is 
invertible. Let U be obtained from U by replacing the odd-positioned vector fac¬ 
tors V r , V r _ 2 ,... counted from the right in the product with J~V r , J r V r _ 2 ,.... 
Let the result be U = W 3 W 2 • • • W r . Then the 4x4 matrix form of Ad\ j is up 
to scale that of the following transformation: 

X£R 4 ^W r (W 2 •(•••• (W r • X))) £ R 4 , (4.15) 

where the inner product is in CZ(R 4 ). 

Definition 14. Let U = X 3 X 2 • • • X r £ Pin{ 3,3). The inverse-transpose of U 
is the following element of Pin( 3,3): 

U J := (JXiXJ-Xa) • • • (JXr). (4.16) 


Lemma 13. Let U £ Pin( 3,3). Then as matrices in 0(3,3), 


(Adu) T — Ad\j a . 


(4.17) 
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Proof. For any M e 0(3,3), M T JM = J. So MT 7 = JMJ. For M = Ad\j, 
for any X e M 3 ’ 3 , 


M~ T X = ./(MIX' 7 )) = Ad^((X J ) J ) = Ad* VJ X. □ 


Below we investigate 3-D projective transformations induced by two reflec¬ 
tions in R 3 ’ 3 . Let U = (xi,yi) T and V = (x 2 ,y 2 ) T , then Ad^j^ has the 
following 4x matrix form up to scale: 


f 0 - Xl T \ / 0 -y 2 T 
\xi yi x I 3x3 / \y 2 x 2 x I 3x3 
xi • y 2 (xi x x 2 ) t 

-yi x y 2 (x 2 • yi)I 3X 3 +xiy 2 


X2yi 



(4.18) 


The resulting matrix has trace 2(x 3 • y 2 + X 2 • yi), so it is in sl( 4) if and only if 

u-v = o. 

When xi = x 2 = x, (4.18) becomes 


( x ■ y 2 0 

V — yi x y 2 (x • yi)I 3x3 + x(y 2 - yi) T 


(4.19) 


For x = ei, and yi 


ziei + z 2 e 2 , and y 2 

= 2/ i e i + 

y 2 e 2 + 2 / 3 e 3 , (4.19) becomes 

/ Vi 

0 

0 

°\ 


-222/3 

2/1 

2/2 - 2 2 

2/3 

, (4.20) 

212/3 

0 

21 

0 

\ 222/1 - 2 i 2/2 

0 

0 

*1/ 



where the submatrix A composed of the last three rows and the last three 
columns is upper-triangular. A is orthogonal if and only if 


j/i= ±2i, y 3 = 0 , 2/2 = 22, 


i.e., either y 3 = y 2 or y 2 = — 2 ie 3 + z 2 e 2 . In the former case, (4.20) equals l 4 X 4 
up to scale; in the latter case, (4.20) becomes 


/ 1 0 0 0 \ 
0 10 0 
0 0-10 
\2z 2 0 0 -1/ 


(4.21) 


The 3x3 submatrix composed of the last three rows and the last three columns 
represents a rotation of angle n about the axis in direction e 3 and passing through 
point z 2 e 3 . Thus we get the following: 

Proposition 22. Let x £ R 3 be nonzero, and y £ R 3 satisfies x • y = 0. The 
rotation of angle it about the axis passing through point y in direction x is gen¬ 
erated by two reflections induced by the following vectors sequentially (from right 
to left): 

x \ / 

Ax + x x y/ ’ (-Ax + xx y 


(4.22) 
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where A ^ 0 and is arbitrary. Except for such rotations in £ 3 , no other rigid-body 
motion is generated by two reflections in R 3,3 . 

There are five kinds of 2-planes in R 3,3 that are spanned by invertible vectors, 
classified by the signature of the 2-plane: R 2,0,0 , R 0,2,0 , R 1 ’ 0,1 , R 0 ’ 1 ’ 1 , R 1 ’ 1 - 0 . We 
use them to give the normal form of any element of <SX(4) that is the product 
of two invertible skew-symmetric matrices, under the matrix similarity transfor¬ 
mations of 5X(4). 

Definition 15. Two matrices M, N £ SL{ 4) are said to be conjugate in SL{ 4), 
if there exists a matrix C £ SL( 4) such that M = CNC' 1 . 

Lemma 14. Let U, U 0 ,V,V 0 £ R 3,3 be invertible, such that U 2 = Uq, and 
U A V has the same signature with Uo A Vo- Then there exists g € SOq(3,3) 
such that g(U) = Uo, and g(V) £ Uo A Vo- 

Proof. Obviously there are two 4-spaces V 4 and Vq of R 3,3 , each having 
signature R 2 ’ 2 , such that U, V G V 4 and Uo, V 0 £ Vq, and the orthogonal map 
h defined on U A V by h( U) = U 0 and h(V) £ U 0 A V 0 can be extended to an 
orthogonal map from V 4 to Vq, still denoted by h. 

Let V 2 be the orthogonal complement of V 4 in R 3,3 , and let V 2 be the or¬ 
thogonal complement of Vq in R 3,3 . Let v + ,v_ be an orthogonal basis of V 2 
such that v^_ = 1 and v 2 = — 1, and let v 0+ , v 0 _ be an orthogonal basis of Vq 
such that Vq + = 1 and Vq_ = — 1. Extend h to an orthogonal transformation of 
R 3 ’ 3 by setting /i(v + ) = vo+ and fo(v_) = v 0 ^. Set g = h on V 4 , and 

— if h is the composition of the reflections with respect to even number of 
positive vectors and even number of negative vectors, set g(v + ) = Vo+ and 
s( v -) = v o-; 

— if h is the composition of the reflections with respect to odd number of 
positive vectors and odd number of negative vectors, set g(v + ) = — vo+ and 
5 (v-) = -v 0 _; 

— if h is the composition of the reflections with respect to even number of 
positive vectors and odd number of negative vectors, set g(v + ) = vo+ and 
3 (v-) = —v 0 _; 

— if h is the composition of the reflections with respect to odd number of 
positive vectors and even number of negative vectors, set g(v + ) = — Vo+ and 
S( v -) = v 0 -. 

Then g is the composition of the reflections with respect to even number of 
positive vectors and even number of negative vectors, so g £ SOq(3,3). EH 

Lemma 15. Let U = X 1 X 2 ...X 2 /C £ Spin( 3,3), where each X,; £ R 3,3 is 
invertible. Let g £ SOo(3,3) be the Pliicker transform of G £ SL{ 4). Denote 
g(U) := (gXi)(gX 2 )... (< 7 X 2 ^). Let the 4x4 matrix forms of Ad{j and ^4d*(u) 
be A, C £ SL( 4) respectively. Then 


A = GCG E 


(4.23) 
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Proof. Denote gX. t = Yj for 1 < i < 2k. Let g = Ad^ for some W £ 
Spino(3,3), then g~ x = Ad^-i- We have 

Ad* v (X) = (g-\ Y 1 Y 2 ■ ■ • Y 2fe ))X(<T 1 (YiY 2 • • • Y 2fe )) -1 

= W(Y : Y 2 • • • Y 2fe )W- 1 XW(Y 1 Y 2 • • • Yafc)" 1 W" 1 

= {Adw o Ad*^^ o Ac?^ v _ 1 )X. |—| 

Proposition 23. Let A be the product of two skew-symmetric matrices of SL( 4). 
Then A is conjugate in SL( 4) to one and only one of the following block diagonal 
matrices up to scale: 


c — s 0 0 

s c 0 0 

0 0 c — s 

0 0 sc 


where c = cos 6, and s = sin# 0. 


(4.24) 


c s 0 0 

—s c 0 0 

0 0 c — s 

0 0 sc 


where c = cos 6, and s = sin# yA 0. 


(4.25) 


1 -A 0 0 
0 10 0 
0 0 10 
0 0 A 1 


where A ^ 0. 


(4.26) 


1 A 0 0 
0 10 0 
0 0 10 
0 0 A 1 


where A ^ 0. 


(4.27) 


diag(A 1 ,A 1 , A, A), where A > 0. 


(4.28) 


diag(A , A 1 , —A, —A), where A > 0. 


(4.29) 


Proof. When A = ±I 4 X 4 , then it belongs to the class (4.28). Below we assume 
A ±I 4 X 4 . Then A must be the 4x4 matrix form of Ad{j V for some unit vectors 
U, V £ R 3 ' 3 , where UAV/0. 

Case 1. If U A V has signature R 2,0,0 , choose 


V n — ■= 


1 ( e 2 


V2 V e 2 


, U 0 = ^ ( 62 + ^ ( 63 ) , (4.30) 
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where c = cos 9, and s = sin 9 ^ 0. They are both unit vectors, and U 0 A V 0 
has signature R 2,0,0 . By Lennna 14, there exists an element g G SO o(3, 3) such 


that g(V) 


V 0 , and g(U) G V 0 A 



Obviously there exists a parameter 9 


such that g(U) = Uo- So the 4x4 matrix form of i s conjugate in SL( 4) 

to that of Ad{j oVo . By (4.18), the 4x4 matrix form of Ad{j oVo is (4.24). 

Case 2. If U A V has signature R 0,2,0 , then choose unit vectors 


V 








(4.31) 


where c = cos 9, and s = sind ^ 0. Obviously U 0 A V 0 has signature R 0,2,0 , so 
an element g G SOq{3, 3) maps V to Vo, and maps U to Uo for some parameter 
9. The 4x4 matrix form of Ad \ JoVo is (4.25). 

Case 3. If U A V has signature R 1,0,1 , then choose unit vectors 


V 




U 




_A_ 

71 



(4.32) 


where A ^ 0. Since U 0 A V 0 has signature R 1,0,1 , an element g G SOq(3,3) maps 
V to Vo, and maps U to Uo for some parameter A. The 4x4 matrix form of 
A^UoVo (4-26). 

Case 4. If U A V has signature R 0,1,1 , then choose unit vectors 


v “ = ^(-7’ u » = 7i(-«) + ^(o)' (4 - 33) 

where A ^ 0. The 4x4 matrix form of 2ldu 0 v 0 (4-27). 

Case 5. If U A V has signature R 1 ’ 1,0 , there are four subcases: (1) U, V are 
both positive, (2) U is negative while V is positive, (3) U is positive while V is 
negative, (4) U, V are both negative. 

Let V' G U A V be a unit vector orthogonal to V. Then V, V' have opposite 
signatures, and U • (V A V') is a unit vector of U A V having opposite signature 
to U. Since 


UV = U(V A V')(V A V')V = -V 2 {U • (V A V')}V', 


we have = Ad*^ .( VAV /)j V /- So we only need to consider the first two 

subcases by assuming V 2 = 1. Then (V') 2 = —1. Let U = aV + bV' where 
b ^ 0. Then |a 2 — b 2 \ = l,soa/ ±b. 

Choose unit vectors 


V 




V(,= 


1 

71 



(4.34) 


An element g G SOq(3, 3) maps V to V 0 , 
up to scale to 


((a + fo)ei\ 
\{a-b)e 1 ) 


and maps V' to Vg, so it maps UV 



(4.35) 
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The 4x4 matrix form of the adjoint action of (4.35) is 


diag(a + b, a + b, a — b, a — b). 


(4.36) 


Let A = 


a—b 

a+6 


. Then A > 0. 


Subcase 5.1. When U 2 = 1, then a 2 — b 2 = 1. The 4x4 matrix form of the 
adjoint action of (4.35) is diag(A _1 , A” 1 , A, A). 

Subcase 5.2. When U 2 = —1, then a 2 — b 2 = —1. The 4x4 matrix form of 


the adjoint action of (4.35) is diag(A 1 , A 1 , —A, —A). 


□ 


Remark. When A ^ 0, then 


(/ix) = (A* + A) + (R ~ A)E(x) A E'(x). (4.37) 

For different parameters A, /i, as long as the ratio A : /i is constant, then 
(4.37) gives the same spinor up to scale. In particular, as long as A ^ 0, then 

^ induces the same rotation of angle 7r about the axis in direction 

x at the origin. 

The following is a general criterion on the “compressibility” of a spinor in 
factored form. 




Proposition 24. Let U be the Clifford product of 4 invertible vectors o/R 3,3 . 
Then U can be written as the Clifford product of 2 invertible vectors of R 3,3 if 
and only if ( U )4 = 0 and (U )2 has an invertible vector factor. 

Proof. Let U = X 1 X 2 X 3 X 4 where each vector Xj is invertible. When U = 
Y 1 Y 2 for some other two invertible vectors Yi, Y 2 , the two conditions of the 
proposition are obviously true. Conversely, when both conditions hold, then the 
X,; are in some 3-space of R 3,3 , so (U )2 must be decomposable into the outer 
product of two vectors. Let a, b be two invertible vector of (U )2 such that 
(U )2 = a A b. We claim that there are two invertible vectors Y t = A^a + /i,b for 
i = 1,2 such that U = YiY 2 . 

Expanding U = (U ) 0 + a A b = (Aia + /zib)(A 2 a + ^b), we get 


A1M2 — A2/J1 — 0 , 

AiA 2 a 2 + /xi^ 2 b 2 + a • b(Ai/i 2 + X 2 R 1 ) ~ (U)o- 


(4.38) 


When setting A 2 = 0, then = A 1 1 and p,\ = ((U ) 0 — a • b)Aib 2 , so U = 
(a + ((U)o — a • b)b _ 1 )b. Since b and U are both invertible, so is Ub” 1 . EH 

In the rest of this section, we construct 4-tuples of reflections in R 3,3 inducing 
rigid-body motions in £ 3 . Affine transformation (4.19) is translation-free if and 
only if yi x y 2 = 0. Consider the composition of two such affine transformations. 
Let 



(4.39) 
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where y 2 = A,z,; for * = 1,2 such that \ ^ For simplicity, assume x,; • z,; = 1 
for * = 1 , 2 . 

The 4x4 matrix form of Ad* is up to scale the following: 

U, ViU 2 v 2 

(o A1I3X3 + (1 ~ Ai)xiZi T ) (o A2I3X3 + (1 — ^2)x2Z2 T ) (0 m) 

(4.40) 

where 


M — AiA2l3x3+Al(l — A 2 )x 2 Z 2 T + A2(1 — Ai)xiZi T + (1 — Ai)(1 — A 2 )(x 2 'Zi)xiZ 2 T . 

(4.41) 

When M is the rotation matrix 

/ cos 6 — sin 6 0 \ 

I sind cos0 0 J , (4.42) 

l 0 01 / 

there is more than one way of constructing the four reflection vectors (4.39). 

Option 1. Choose xi = ei and X 2 = e 2 - Let zi = ei + 0 , 2^2 + 0363 and 
Z 2 = 6 iei + e 2 + & 3 e 3 . Comparing (4.41) and (4.42), we get 


Ai = cos 9, X 2 = sec 9, 



So 


Ux 

Vx 

U 2 

V 2 


ei 

— e 9 cos 9 cot | 

• t.\ " 

cot I / 


e 2 

- ei sec 9 cot f 


(4.43) 


(4.44) 


This 4-tuple of reflection vectors has the strong defect that the following is 
required: 9 ^ {0,±7r/2}. 

Option 2. Set xi = ei and X 2 = c'e-i + s'e 2 , where c ' 2 + s' 2 = 1. Set 
zi = ei + ae 2 and z 2 = pei + qe 2 such that x 2 • z 2 = cp + sq = 1. Set yi = Azi 
and y 2 = A _1 Z 2 for some A / 1. Again by comparing (4.41) and (4.42) and 
solving the equations, we get 


A = -1, 


a = 0 , 


s =q = 


— sin—, c =p = cos-. 
2 2 


(4.45) 
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So 


ei 

-ei 

ei 

ei 


ei cos | — e 2 sin | 
—ei cos | + e 2 sin | 
er cos I — e 2 sin | 


(4.46) 


Ur = 

V-, = 

U 2 = 

Vq = . 

ei cos 2 — e 2 sin ^ 

This 4-tuple of reflection vectors is perfect. It can be reformulated as follows: 


Proposition 25. Let xi,x 2 be unit vectors o/R 3 . The rotation of R 3 in the 
xix 2 plane with angle 2Z(xi,x 2 ) is generated by four reflections in R 3,3 with 
respect to the following invertible vectors sequentially (from right to left): 


x 2 

-x 2 



(4.47) 


Remark. The adjoint action of 




induces the rotation of angle 7r 


in the plane X] 1 . It is a geometric fact that this rotation followed by the rotation 
in the plane x^ leads to the rotation of angle 2 Z(xi,x 2 ) in the Xix 2 plane. 

Next consider the translation by vector t. In the 4 vectors of (4.39), denote 
w » = y i — z i for i = 1,2, and assume x* • z$ = 1, then the matrix form of 
is up to scale the following: 


1 0 
-Wi X Zi (Xi • Wi + 1 ) 13 x 3 - XiWi T 
1 0 
fi h 

where 


1 0 

-w 2 X Z 2 (x 2 • W 2 + 1 ) 13 x 3 — x 2 w 2 t 


(4.48) 


/l = — wi x zi — (xi • wi)w 2 x z 2 — w 2 x z 2 + xi[wiw 2 z 2 ], 
fi = (xi • Wi + l)(x 2 • w 2 + l)l 3x 3 - (xi • w i + 1)x 2 w 2 t (4.49) 

— (x 2 • W 2 + l)XlWl T + (x 2 • Wl)xiW 2 T . 


Notice that in (4.49), vectors zi,z 2 occur only in the translation part f\. 
That (4.48) is a translation matrix if and only if / 2 = 13 x 3 - For simplicity, 
choose xi = x 2 = x, and set x = ei. 

Option 1. To simplify the expression of / 2 , let x ■ w, = 0 for i = 1,2. Then 
equation / 2 = I 3 X 3 becomes 

—x(w 2 + (x • wi)w 2 — Wi) T = 0, (4.50) 


hence wy, w 2 are linearly dependent, and we get w 2 = — wy. Choose wi = e 2 . 
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The translation vector is t = /i = — Wj x (zi — Z 2 ). We can choose t to be 
in the direction of ei, so that t = — Aei for some scalar A. Choose Z 2 = ei, then 
zi = ei + Ae 3 . The 4 reflection vectors inducing the translation by vector — Aei 
are the following from right to left: 


ei 

ei + Ae 3 


e 2 


ei 

ei + Ae 3 


ei 

e 3 — e 2 



(4.51) 


This 4-tuple can be formulated in a coordinate-free form as following: 

Proposition 26. Let x, y be nonzero vectors of R 3 satisfying x • y = 0. The 
translation of R 3 by vector x is generated by four reflections induced by the 
following invertible vectors in R 3,3 sequentially (from right to left): 


x- y + 


X x y 


x 

X X y 


x 

x + y 


(4.52) 


Remark. The first two reflections generate the following affine transformation: 

1 0 
-xxy + x I 3x 3 + xy T 

the last two reflections generate the following affine transformation: 

1 0 
xx y I 3x3 -xy T 


(4.53) 


(4.54) 


The two matrices are always commutative. 

Geometrically, when setting |y| = 1, then (4.54) represents a planar shear 
transformation followed by a translation by vector x x y, where the shearing 
occurs in plane (x x y) 3 -, with shear direction x and shear factor — |x|: for any 
z £ R 3 , z 1—^ z (z • y)x + x x y. Similarly, (4.53) represents a planar shear 
transformation followed by a translation by vector x — x x y, where the shearing 
occurs in direction x of plane (x x y) 3- , while the shear factor is |x|. 

Option 2. Choosing x ■ w, = —2 for i = 1,2 simplifies the expression of / 2 
just the same. Then equation / 2 = I 3X 3 becomes 

—x(w 2 + (x • w 3 )w 2 + w 3 ) t = 0, (4.55) 


hence wi,w 2 are linearly dependent, and we get w 2 = wq. Choose wq = w 2 = 
—2 e 3 . 

The translation vector is t = f\ = —wq x (z 3 — z 2 ). We can choose t to 
be in the direction of e 3 , so that t = —Ae 3 for some scalar A. Choose z 2 = e l5 
then zi = e 3 — Ae 2 /2. The 4 reflection vectors inducing the translation by vector 
—Ae 3 are the following from right to left: 


ei 

A 

—e 3 — —e 2 


ei 

A 

e 3 — —e 2 




(4.56) 


This 4-tuple is just as perfect as the previous option. Its coordinate-free formu¬ 
lation is the following: 
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Proposition 27. Let x, y be nonzero vectors of R 3 satisfying x • y = 0. The 
translation of R 3 by vector x is generated by four reflections induced by the 
following invertible vectors in R 3,3 sequentially (from right to left): 


x x y 
2 


(4.57) 


Remark. The first two reflections generate the following affine transformation: 


(4.58) 


x — I 3 x3 + 2- 


it is the rotation of angle 7r about the axis in direction y at point x/2 £ 
last two reflections generate the following affine transformation: 


0 -I 3 


It is the rotation of angle 7r about the axis in direction y at the origin. The two 
matrices are not commutative. 

We work out the spinor representation in factored form of a general rigid- 
body motion (R, t) : x 6 R 3 i— >• Rx + t, where R is a rotation matrix, and 
t € R 3 . It is a classical result that a rigid-body motion in space is either a pure 
translation or a screw motion: a rotation about an affine line called the screw 
axis, and then a translation along the screw axis by a screw driving distance. 

Lemma 16. Let R be a rotation of angle 6 about the axis in unit direction v 
at the origin, and let t • v = 0. Then the rigid-body motion x i—> Rx + t of R 3 
is a pure rotation, and the axis of rotation is the line in direction v and passing 
through point 

Rtt- 9 t 

c=—(4-60) 
2 sin | 

where R Q denotes the rotation of angle a about unit direction v at the origin. 

Proof. We only need to prove that c is a fixed point. Set 

( cos 9 sin 9 0 \ 

— sin# cos 9 0 . (4.61) 

0 0 1 / 


Substituting them into c — Rc = t, we get 



where c = cos(0/2) and s = sin(0/2). The solution is 



(4.62) 

(4.63) 
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Corollary 4. Let R be a rotation of angle 9 about the axis in unit direction v 
at the origin. Then the rigid-body motion x i —> Rx + t of R 3 is a screw motion, 
whose screw axis is in direction v and passing through point 


R, tt — 6 (t — (t • v)v) 
2 sin | 


(4.64) 


In fact, point c is the foot drawn from the origin to the screw axis, called the 
original center of the screw motion. The screw driving distance is d = t • v. 


Proof. Since t = (t — (t • v)v) + (t • v)v, point c of (4.64) is first mapped to 
Rc + (t — (t • v)v) = c by (4.60), then mapped to c + (t • v)v, so it is on the 

screw axis. EH 


Consider a rigid-body motion (R, t) where the rotation axis of R is e 3 . Let 
ei be rotated to ei cosd + e 2 sin 9, and let t = t\&i + t 2 e 2 + £363 be the transla¬ 
tion vector. Motivated by the 4-reflection generators of pure rotation and pure 
translation, consider the 4 reflections generated by the following vectors: 


Ur = 


where 


Vr = 


U 2 = 


cei — se 2 

qi 


v 2 - 


cei — se 2 

q2 


(4.65) 


(4.66) 


Yi = ~ e i + o: 3 e 2 — a 2 e 3 , 
y 2 = ei + a 3 e 2 — a 2 e 3 , 
qi = -(cei - se 2 ) + /3 3 (se 1 + ce 2 ) - /3 2 e 3 , 
q 2 = (cei - se 2 )+/3 3 (sei + ce 2 ) -/I 2 e 3 . 

The 4x4 matrix of the affine transformation induced by the above 4 reflections 
is 


( 1 0 H 

V —2(a 3 e 3 + a 2 e 2 ) —I 3x3 + 2eief y 

( 1 0 H 

V — 2(^ 3 e 3 + /? 2 (se i + ce 2 )) — I 3X 3 + 2(cei — se 2 )(cei — se 2 ) T J 


1 0 
M R' ' ’ 


where 


0 


By t' = t, we get 


cii + st-2 ti t 3 

«2 —-T-, p 2=— — , -Ol3+P3 = —. 

Zs Zs Z 


(4.67) 


t' = ei(—2 s/3 2 ) + e 2 (—2a 2 + 2 c/3 2 ) + e 3 (—2a 3 + 2/3 3 ), 

' cos 9 — sin 9 O' 
R' = I 3 X3 — 2eief + 2(cei + se 2 )(cei — se 2 ) T = I sind cos 9 0 


0 1 


(4.68) 

(4.69) 
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When s/0, i.e., 9 0, we can set 


0 3 = 0 , 0 3 = — 


h 


(4.70) 


Then by (4.64) and d = < 3 , 


d 


= e i x (c+ -e 3 ), 


0:362 — 02e 3 
/3 3 (sei + ce 2 ) - /3 2 e 3 = (ce 3 - se 2 ) x c. 

Alternatively, we can set 


Then 


0 3 e2 — 02 e 3 


03 = 0, 03 = . 


= e 3 x c. 


/S 3 (sei + ce 2 ) - 0 2 e 3 = (ce 1 - se 2 ) x (c - -e 3 ). 


(4.71) 


(4.72) 


(4.73) 


When s = 0, the rigid-body motion is a translation, and there is no rotation 
plane eie 2 . Still as long as t\ = 0, we get from (4.69) and (4.72) /3 2 = o 3 = 0, and 
02 = —£ 2 / 2 , and 03 = t 3 /2. The 4 reflection vectors generating the translation 
are (from right to left): 


ei 


-er 


h 


ei 


ei 


-e 3 


ei 


-e 3 


-ei 


* 3 


-e 2 


ei 

i 3 

, ei + y e! . 


(4.74) 


It gives the third 4-tuple of reflection vectors generating pure translation by 
t = t 2 e 2 + t 3 e 3 other than (4.51) and (4.56). 

Proposition 28. Any rigid-body motion is induced by two pairs of reflections 
in K 3,3 such that each pair induces an affine transformation. For a rotation in 
a plane spanned by two orthonormal directions vi, v 2 with angle 6 0, followed 

by a translation along vector t, the corresponding 4 reflections are the following 
(from right to left): 


cvi + sv 2 
— (cvi + sv 2 ) + mi 


CVi + sv 2 
cvi + sv 2 + mi 


Vl 

—vi + m 2 


where c = cos( 0 / 2 ) and s = sin( 0 / 2 ), and we can choose either 


Vl 

Vl + m 2 

(4.75) 


mi = (cvi + sv 2 ) x (c + -n), 


m 2 = Vi x c, 


or 


mi = (cv 1 + sv 2 ) x c, m 2 = Vi x (c - -n). 


(4.76) 


(4.77) 


Here c is the original center (4-60) of the screw motion, and d is the screw 
driving distance. 
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By (4.76), the first two reflections generate the rotation of angle tt about the 
axis in direction cvi + sv 2 at point c + dn/2, and the last two reflections generate 
the rotation of angle i r about the axis in direction vi at point c. By (4.77), the 
first two reflections generate the rotation of angle n about the axis in direction 
cvi + sv 2 at point c, and the last two reflections generate the rotation of angle 
7 r about the axis in direction Vi at point c — dn/2. When 0 = 0, for the choice 
of vi such that t • Vi = 0, (4.75) is still valid by taking the limit of 0 —> 0. 

Corollary 5. Any rigid-body motion is the composition of two 3-D rotations 
of angle 7r. 

Given a 3-D projective transformation or polarity, we see that the elements 
of Pin sp ( 3,3) generating it in factored form is far from being unique. On the 
other hand, the expanded forms of the elements differ at most by a factor in 
{±l,±I 3j3 }. For a rigid-body motion, the exponential of a bivector of /1 2 (R 3 ’ 3 ) 
provides a spinor that when expanded is equal to the spinors in different factored 
forms that we find in this section. This is the topic of the bivector representation 
of se(3), the Lie algebra of rigid-body motions. 


5 Lie subalgebra of rigid-body motions and classical 
screw theory 


First we recall some basic facts of the Lie algebra se(3). By direct verification, 
we know that 


Proposition 29. Let R be the rotation matrix of angle 0 about the axis in unit 
direction v at the origin. Then 


(1 r)=“»(u OvxVj’ 


(5.1) 


where 


u = (t • v)v + 


2 sin 


fR_ e (t — (t • v)v) = 0C X V + dv. 


(5.2) 


and the original center c is given by and d = t • v. 

The Lie algebra se(3) is composed of 4x4 matrices of the form f ^ 

v x I 3x3 

where t, v € R 3 . Hence we can use a pair of vectors t, v to represent the matrix. 


Definition 16. The map 

(2 vxl 3x3 ) eSe(3W (t) GR3 ’ 3 (5 ' 3) 

is called the screw representation of se(3). The Lie bracket of se(3) agrees with 
the following cross product of screw forms: 

f Vl ) X f V2 ) ■= ( V 1 X V 2 

VW V* 2 / V V 1 X t 2 + V 2 x ti 


(5.4) 
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By (5.1), the screw form of the element of se(3) generating the motion x £ 
R 3 i—i Rx +1 where R ^ 13 x 3 , is 


0 ( " ), (5.5) 

(cxv + Av) v ’ 

where A = d/6 = (t • v)/9 is called the screw ratio. It is the ratio of the screw 
driving distance d by the angle of rotation 6. The screw motion is positive if 
A > 0, in which case the screw driving direction and the rotation orientation 
follow the right-hand rule. The screw motion is negative if A < 0, in which case 
the screw driving direction and the rotation orientation follow the left-hand rule. 
The screw motion is a pure rotation if A = 0. 

Any nonzero vector X £ R 3,3 is up to scale either of the form 

(t)=E'(t), (5.6) 

or of the form 

X =(cxl‘+Al)' < 5 ‘> 

where 1 is a unit vector of R 3 , c £ R 3 and /i € R. In terms of infinitesimal 
generators of screw motions, (5.6) generates a pure rotation in direction t. As 
to (5.7), since X 2 = 2A, the vector is null if A = 0, and represents a line in 
direction 1 and passing through point c, at the same time it is the generator of a 
pure rotation about the axis represented by the 6-D vector itself. The vector is 
positive if A > 0, and generates a positive screw motion about the axis ( 1 , c x 1 ) T ; 
the vector is negative if A > 0, and represents a negative screw motion about 
the same axis. So the screw representation of se(3) provides a motion-generator 
interpretation to all vectors of R 3,3 . In contrast, in line geometry only null vectors 
of R 3,3 have geometric interpretation. 

Traditionally, a “screw” refers to the screw form of a pure rotation or pure 
translation, i.e., a null vector of R 3,3 , while a general vector of R 3,3 is called a 
twist, as it generates a screw motion. In this paper, we call the vectors unani¬ 
mously a screw form. 

The following proposition states that the cross product of the screw forms of 
two screw motions is the screw form of a third screw motion, the latter being 
along the common perpendicular of the two screw axes of the former two motions. 


Proposition 30. Fori = 1,2, let 


X, = 


h 

c / x 1, T A, 1 j 


i>3,3 


(5.8) 


such that (1) each 1* is a unit vector, (2) c • 1; = 0, and (3) li x b ^ 0. Let 
y £ R 3 he an arbitrary point on the common perpendicular of the two lines 
(lj, Cj x lj), and let d be the signed distance from the first line to the second along 
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the direction li • I 2 , then 


X, x X 2 = 


where 


li x I 2 

y x (li x I 2 ) + /-tli x I 2 

I1I2 


(5.9) 


/i — Ai + A 2 + d- 


li x 1 2 | 


(5.10) 


For example, the following point of R 3 is on the common perpendicular: 

(Cl ' 12)12 + (C2 ’ ll)ll 


y = 


and 


d = 


(li x 1 2 ) 2 

[(c 2 - CiJUa] 


111 X 1 2 | 

Proof. Obviously, li, I 2 , li x O form a basis of R 3 . Let 

c i = P 2 (I 2 — (li • I 2 ) 1 1 ) + Pill x I 2 , 
c 2 = <72(ll ~ (ll ' 12)12) + 91^1 X I 2 . 

Then (5.11) and (5.12) become 

y = P 2 I 2 + 92 I 1 , d = (qi - pi)\h x l 2 \ 


(5.11) 

(5.12) 

(5.13) 


It is easy to verify that line 


meets both lines 


1 , 


C; x 1, 


(5.14) 
for 


li x 1 2 

y x (h x 1 2 ) 

* = 1 , 2 , as the inner products of the vector with the other two are both zero; so 
the line is the common perpendicular of the other two lines, and d is the signed 
distance between the other two lines. 

By direction computation, we get that Xi x X 2 equals 

li x 1 2 \ 

li x l 2 (Ai + A 2 + (<71 — Pi)li • 1 2 ) + li(_P 2 + 92 I 1 • I 2 ) — 12(92 +P 2 I 1 • I 2 )) 


By y x (li x 1 2 ) = li(p 2 + 92 I 1 • I 2 ) - 12(92 + P 2 I 1 • I 2 ), we get (5.9). □ 

Proposition 31. 

((qi) X (q 2 )) ' (q 3 ) = t 11 * 2 ^ + [ l2l3<l1 ^ + t l3liq2 ]' ( 5 ' 15 ) 

The result is invariant under a shift of the subscripts 1,2,3. In particular, when 
q, = c i x I, + Ajlj, the result is 


ll' (C 3 — C 2 XI 2 - 13 ) + I 2 • (Ci — C 3 )(l 3 - li) + I 3 • (C 2 —Ci)(li T 2 ) + (Ai -|- A 2 + A 3 )[I 1 I 2 I 3 ]- 

(5.16) 
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Lemma 17 (Decomposition of vectors of R 3,3 with respect to the cross 
product). For any x,y G R 3 , let z G R 3 be nonzero such that z • x = 0 and 
z • y = 0. Then 




(5.17) 


Proposition 32 (Orthogonal decomposition of infinitesimal screw mo¬ 
tions of R 4 ). Let v i,V 2 ,V 3 be an orthonormal basis of R 3 . Then 


Vl ^ = ( V2 ^ x ( V3 
Av 3 x Vl + /XV1 y \ Av 3 x v 2 y \ /./.V 3 

/t states that an infinitesimal screw motion about the axis in direction vi and 
through point Av 3 , is the cross product of an infinitesimal rotation about the axis 
in direction v 3 x v 3 and through point Av 3; with an infinitesimal screw motion 
about the axis in direction v 3 and through the origin. 

We come back to the definition of the cross product in R 3,3 . It has some 
limited covariance in 50(3,3). Consider B € 50(3, 3). Let its matrix form with 

respect to Witt decomposition R 3,3 = I 3 ® J 3 be ( ^ 12 ). Matrix B is called 

\b 2 i b 22 y 

affine if its induced transformation in R 4 preserves the plane at infinity of 5 3 , 
i.e., B 32 = 0. When B is affine, by B T JB = J, we have 

B 22 = B n T , B^Bu = -B^B 21 . (5.19) 

When B is affine, it is Euclidean if its induced transformation in R 4 preserves the 
inner product of R 3 = (ei, e 2 , e 3 ), i.e., matrix Bn is orthogonal, and B 22 = Bn. 
If B induces a translation, then Bn = B 22 = ±I 3x3 ; if it induces a rotation 
about the origin, then B 23 = 0. 

Proposition 33. The cross product is covariant under the Euclidean subgroup 
of 50(3,3) inducing rigid-body motions of £ 3 , i.e., for any 

M= (c r) €50(3 ’ 3) ’ (5 ' 20) 


(5.18) 


where R is a rotation matrix, and C = —AC 1 A, 


Proof. By direct computation, we get that the cross product is covariant under 
the subgroup of 50(3, 3) inducing translations, and is also covariant under the 
subgroup of 50(3,3) inducing rotations at the origin. So it is covariant under 

the Euclidean subgroup. EH 

In classical screw theory, a force is represented by a line in space, and a torque 
is represented by a line at infinity. For example, the force along line (ep + x)l 
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with magnitude / is represented by /(e 0 l + xl) £ T 2 (R 4 ), the torque in the plane 
normal to unit vector n with scale t, is represented by £ /1 2 (R 4 ). The linear 
space spanned by forces and torques is the space of wrenches. A general wrench 
is represented by a vector of R 3,3 , while a pure force or torque is represented by 
a null vector of R 3 ’ 3 . 

When an infinitesimal rigid-body motion is represented in screw form, i.e., 
as a vector of R 3,3 , then the virtual work of a wrench and a screw form is defined 
as the inner product of the two corresponding vectors in R 3,3 . Let (f, q) T be a 
wrench, where p is the force part and q is the torque part, and let (v, u ) T be 
the screw form of a screw motion with axis direction v. Then the virtual work 
between them is 




f ■ u + q • v. 


(5.22) 


We try to under the virtual work in the 4x4 matrix form of se(3). By 
(4.12), any vector X £ R 3,3 corresponds to a bilinear form £(X) £ sl(4). When 
X = (f,q) T represents a wrench, then the inner product of £(X) with the 

infinitesimal motion matrix (° ° T ) is 

Vu vxl 3x3 / 


(° - fT 'WO 0 \ 

V f q X 13 x 3 / V U v X 13 x 3 / 


—f • u — 2(q • v). 


(5.23) 


It is almost the virtual work except for the ratio of the two part of works by the 
force and the torque respectively. 

Since Spin( 3,3) covers SL(4), the Lie algebra A 2 (R 3,3 ) of Spin( 3,3) is iso¬ 
morphic to si(4). In particular, any element of se(3) has a bivector form in 
T 2 (R 3 ' 3 ). On the other hand, a wrench is only a vector of R 3,3 . How is the 
virtual work represented by some pairing between a bivector and a vector of 
yl(R 3,3 )? Before embarking on the investigation of this problem, we make an 
analysis of the bivector representation of se(3) by computing the expanded form 
of the spinors (4.46), (4.51), (4.75) inducing rigid-body motions and the corre¬ 
sponding exponential form. 


Proposition 34. Let x = ei and y = cei + se 2 , where c = cos(0/2) and 
s = sin(0/2), then 

K4)©(- x )(x)-^ e -- e -»- <»■”> 

Proof. On one hand, 



= ((cEi + sE 2 ) A (cE; + sE , 2 ))(E 1 A Ei) 
= c 2 + cs(E 2 i' — Ei 2 /) — s 2 Ei 2 i' 2 '. 


(5.25) 
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On the other hand, let B = (E 2 p — Ei 2 ')/2, let I = Ei 2 / 2 i/, and let L = (l+I)/2. 
Then 

I 2 =1, 

BI = IB = B, 

B 2 = -L, 

L 2 = L. 


(5.26) 


So 


exp($B) = 1 + B 9 + 


—L0 2 -B6» 3 L 6 4 Bd 5 


3! 


4! 


5! 




2 ! 

n3 fjb n2 a4 

= B(0 -¥ + 5! + - ) + 1 - L + L(1 -¥ + 4f + - ) 
= B sin 9 -\ ---1-4— cos 9 


2 2 
9 
2 


2 9 . 2 @ 0 . 9 

= cos — — I sin — + 2B cos — sin —. 


□ 


Proposition 35. The translation along vector x £ R 3 is induced by the follow¬ 
ing spinor ofCl(3,3): 


1 _ E(x) J 3 = e xp(_l E (x) • J 3 ). (5.27) 

Proof. Expanding the Clifford product of the four vectors in (4.51) from left 
to right, we get 4 + 2AE 2 / 3 /. The exponential form of spinor 1 + AE 2 / 3 //2 is 

obviously exp(AE 2 / 3 //2). [H 

Definition 17. Let x = (xi, x 2 , x 3 ) £ R 3 . Define the following linear maps from 
R 3 to yl 2 (R 3 - 3 ): 


EE : x = xiei + x 2 e 2 + x 3 e 3 H> x 3 E 23 + x 2 E 3i + x 3 E i2 ; 

E'E' : x = xiei + x 2 e 2 +x 3 e 3 h> XiE 2 / 3 /+ x 2 E 3 /i/+ x 3 Ei/ 2 /; 

E<E' : x = Xiei + x 2 e 2 + x 3 e 3 i —> xiE 23 / + x 2 E 3 i/ + x 3 Ei 2 /; 

E>E' : x = xiei + x 2 e 2 + x 3 e 3 H> XiE 32 / + x 2 E i3 / + x 3 E 2i /. 

Obviously, EE(x) = E'(x) • I 3 , and E'E'(x) = E(x) • J 3 . Under a special 
orthogonal transformation of R 3 , the basis ei,e 2 ,e 3 is changed into another 
basis, and the basis vectors E,, E' of R 3,3 change accordingly. With respect to 
the new basis, four new maps can be defined as above. The following invariance 
can be proved easily: 

Proposition 36. The maps EE, E'E' andEiyE'— are all invariant under 

SO( 3). Furthermore, for any xi,x 2 G R 3 , 


(E>E' - E<E')(x 1 ) • (E>E' - E<E')(x 2 ) = -2( Xl • x 2 ). 


(5.29) 
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The elements of E'E'(x) for x G R 3 form a Lie subalgebra of sl( 4); it is the 
Lie subalgebra R 3 of translations. The elements of (E>E' —E<E')(x) for x G R 3 
form the Lie subalgebra so(3). It is easy to prove that any nonzero element of 
so(3) in bivector form is not the outer product of two vectors of R 3,3 . 

Consider the expanded form of the Clifford product of the 4 vectors in (4.75): 


Eli' —^-Ei/ 2 ' 


c't\ — s'f 2 
2s' 


E 


3'r 


c' Ehi + s 2 E 2 2' + c , s , (Ei 2 / + E 2 p) 

c'ti ti 

-\ -E3/1/ — —Eo's' 

2s' 2 23 


— 4 (c 1 — s 2 E 12 / 21 / — -{sti + ct 2 )En/ 2 / 3 / — -(sf 2 — cti)E 2 i/ 2 / 3 / 

+cs(Ei 2 / — E 2 p) —^-Ei ' 2 ' + —( st 2 — cfi)E 2 / 3 / — —(sti + ct 2 )E 3 /i/^ 

(5.30) 

After some computation, we get that (5.30) equals exp(B 2 ), where if 9 ^ 0, 
then 


B 2 --(E 2 p — Ei 2 /) — -^-Ei/ 2 / 


4 sin ■ 


6 9 9 

(fi cos | + t 2 sin -)E 2 / 3 / + (- ti sin - + t 2 cos -)E 3 / r 


(5.31) 

If 9 = 0, the above expression is still valid by (5.27), as the limit of the right 
side when 6 —> 0 is — (tiE 2 ' 3 / + f 2 E 3 /p)/2, and ei,e 2 are selected so that the 
translation vector t is in the eie 2 plane. 


Proposition 37. The rigid-body motion x 1 —> Rx +1, where R is the rotation 
matrix of angle 6 about the axis in unit direction v at the origin, and t G K 3 , 
corresponds to the following spinor o/R 3,3 , where c is the original center defined 
by ( 4 . 64 ), and d = t • v is the screw driving distance: 

exp i {(E>E' - E<E / )(0v) - E'E'(0c x v + dv)} . (5.32) 

The Lie algebra se( 3) in bivector form is spanned by E'E'(x) and (E>E' — 
E<E')(y) forx, y G R 3 . 

Proposition 38. The cross product o/R 3,3 is induced by the following trivector 
of A 3 { R 3 ’ 3 ): 

C 3 = .F(ei) + .F(e 2 ) + .F(e 3 ) = E 12 / 3 / + E 23 /p + E 31 / 2 /, (5.33) 

such that for any X, Y G R 3,3 , 

X x Y = -(XA Y) • C 3 . (5.34) 

Furthermore, for X = (x,y) T , 

X • C 3 = E'E'(y) - (E>E' - E<E')(x) G se(3). 


(5.35) 
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Let L 2 = (E>E' — E<E')(x) — E'E'(y) be an element of se( 3). The map 



(5.36) 


gives the screw form of bivector L 2 £ se(3). It is also called the se(3)-lift. of 
vector x, y) T . 

The transformation X 1 —L 2 • X for X £ R 3,3 has the following 6x6 matrix 
form: 


M := 


/x x I 3x3 0 \ 

\y x ! 3 x 3 xxi 3x3 y 


(5.37) 


For X = ( P | £ R 3 ’ 3 , 

\q ) 


L 2 • X = MX = f XXp ) = x f P V (5-38) 

Vxxq + yxpy \yj \qj 

So the cross product of two vectors of R 3,3 equals the inner product of the 
se(3)-lift of the first vector with the second vector. When both vectors are lifted 
to bivectors by the se(3)-lift, then the cross product of the two bivectors of se(3) 
equals the cross product of the two original vectors of R 3,3 . 

Proposition 39. By the following image of C 3 under A 3 J: 

D 3 := F"{e\) + P"(e* 2 ) + P"(e* 3 ) = E r23 + E 2 ^ 3 i + E 3 , 12) (5.39) 

we have 

X • D 3 = EE(x) - (E>E' - E<E')(y) £ so(3, 0,1). (5.40) 

Furthermore, 

((u)' C ») ((q)' Ds ) = - f ' U - 2(qV) ' (541) 

Proof. That X ■ D 3 £ so(3, 0,1), the Lie algebra of the special orthogonal 
group of R 3,0,1 , will be made clear in the next section. All others are by simple 

computation. EH 

It is easy to prove that both C 3 and D 3 are invariant under any transforma¬ 
tion of 50(13) and the induced transformation of SO{ J 3 ) with the same matrix 
form. In this sense, the two trivectors are basis-independent. 

By (5.23) and (5.41), for two vectors X, Y £ R 3,3 , if X is lifted by C 3 to a 
bivector A 2 of se(3), and the other is lifted by D 3 to a bivector B 2 of so(3,0,1), 
then their inner product in A 2 (R 3,3 ) is the inner product of the 4x4 matrix 
form of A 2 £ sZ(4) with the bilinear form £(Y). The lift by D 3 realizes a pairing 
between a vector and a bivector of se(3), and the result is almost the virtual 
work between a wrench and an infinitesimal screw motion. 

We have the following more general result on the pairing: 
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Proposition 40. For C 3 ,D 3 given by (5.33) and (5.39), 


((y:)-< C ’ + AW )(ft) (Dl + AJs) ) 


—2xi • y 2 - (1 + A 2 )yi • x 2 . 

(5.42) 


Having talked enough of se( 3), we consider the whole group of Euclidean 
transformations, i.e., affine transformations of the form x £ R 3 i —> Gx + t, 
where G £ 0(3,3). For matrix 


we have 


T 



T:=diag(-1,1,1,1), 


(5.43) 



(5.44) 


Both results are in 5*0(3, 3), so we can simply compute the spinor generator of 
each matrix, say Up U r respectively, and then get the generator ofxi-i Gx +1 
in TSpin( 3,3): TUi = U r T. 

For example, consider the 3-D reflections with respect to affine plane (n,d), 
the plane normal to unit vector n and whose distance from the origin along n is 
—d. The foot drawn from the origin to the plane is f = —dn. Let y £ R 3 . Then 
the image of the reflection of y with respect to the plane is y — 2((y — f) • n)n = 
y — 2(y • n)n — 2dn. The 4x4 matrix form of this transformation is 

M =(-L < 5 - 45 > 


Lemma 18. The screw motion of rotation angle tt and screw driving distance 
d about the axis in normal direction n and passing though point cel 3 , where 
c ■ n = 0, has the following 4x4 matrix form: 


( 1 

y 2c + dn 


0 

2nn T - I 3X 3 


(5.46) 


Proof. By (4.64), the translation vector is t = 2c + dn. The rotation matrix 
of angle tt about the axis in direction n at the origin is 2 nn T — I 3X 3 - EH 

By (5.44), (5.45) and (5.46), we get 

Corollary 6. The reflection with respect to plane (n,d) is the composition of 
T with the screw motion of angle tt and screw driving distance 2d about the axis 
in direction n at the origin. 

Proposition 41. The screw form of se(3) realizing the reflection with respect 
to plane (n, d) is (in, 2dn) T £ R 3,3 . The corresponding element of TSpin(3,3) 
is the following: 


T 


v 2 \ ( V 2 \ / Vl \ / Vl \ 
-V 2 dvi J y V 2 dv 1 J y Vi J \ V 3 J 


(5.47) 
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where vi, V 2 , n form an orthonormal basis o/R 3 . When d = 0, the above element 
can be simplified to 


T 


n 

—n 


(5.48) 


6 Other Lie subalgebras and corresponding screw forms 


The 6-D Lie algebra se(3) of rigid-body motions provides an interpretation of 
vectors of R 3,3 as screw forms of infinitesimal rigid-body motions. The inner 
product of R 3,3 can be transferred to se( 3), although incompatible with the inner 
product of the bivector representation of se(3); conversely, the cross product of 
se(3) in bivector form can be transferred to R 3,3 , and the result is the cross 
product in classical screw theory for line geometry. The transformation-generator 
interpretation of vectors of R 3,3 is not necessarily restricted to the specific Lie 
subalgebra se(3). As long as there is a 6-D Lie subalgebra of sl( 4), a screw 
representation may be assigned to the Lie subalgebra, and the Lie bracket of the 
Lie subalgebra may be translated to a new cross product of the screw forms. 

First we investigate several typical 3-D Lie algebras of sZ(4). 

1. Perspectivity group: 

A perspectivity [13] has the following matrix form: 


(1 x T 

^0 I 3x3 


( 6 . 1 ) 


It induces the following rational linear transformation on the line through the 
origin in unit direction 1: 

A1 € R 3 ^ - h. -1 g R3_ (6.2) 

1 + Ax • 1 

When A —> oo, the point at infinity of the line is mapped to point l/(x • 1). When 
x • 1 = 0, every point on the line is invariant. The origin is fixed, so is every line 
through the origin. 

Every point on the plane normal to x and through the origin is fixed. Fur¬ 
thermore, draw a family of parallel planes normal to x, and let A be the signed 
distance of the plane from the origin along x. Then the plane with signed dis¬ 
tance A is mapped to the plane with signed distance A/(l + A d). In particular, 
the plane with signed distance —1/d is mapped to the plane at infinity, while 
the plane at infinity is mapped to the plane with signed distance 1/d. 

Proposition 42. The bivector generator of (6.1) in A 2 (R 3,3 ) is E'E'(x)/2. The 
Lie algebra of the group is isomorphic to R 3 . 

2. Anisotropic dilation group: 

The transformation (1, ti, t 2 , ^ 3 ) £ GL( 4) where t^t^ > 0, is called an 
anisotropic dilation (or non-uniform scaling). When all the tj > 0, the dila¬ 
tion is said to be positive. All anisotropic dilations of £ 3 form an Lie group, 
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denoted by AD{3)\ all positive anisotropic dilations of £ 3 form a Lie subgroup, 
denoted by AD + (3). We have 

AD(3) = AD+(3) U diag(l, 1, — 1, — l)AD + (3) 

U diag(l,-l,l,-l)AD+(3) U diag(l, — 1, — 1, l)AD + (3). 

Lemma 19. Bivector — AEm/2 induces the following transformation of £ 3 : 

diag(l, 1, e A , e A ). (6.4) 

The corresponding element of si( 4) is diag(—f, j)- 

Proof. Let U = exp(^Em) = cosh ^ + Eu< sinh Then Ad{jEi = e A Ei and 
Ad{ jEp = e- A Ep, while all other basis vectors are fixed. By direct computing, 
we get that Ad\j has the following 4x4 matrix form: diag(e^, e ?, e ~^, e - ^). 

□ 


Definition 18. Denote 

E22' + E33/ — Eh/ tt, E33' + Eh' — E22' _ Eh' + E22' — E33' 

Fi = -. F 2 =-. F 3 =-. 

2 2 2 

(6.5) 

The following are linear maps from R 3 to A 2 (R 3,3 ): 

E = E' : (x±, X2, X3 ) t 1 ^ a:i(E22' + E33') + :C2(E33' + Em) + ^(Em + E22'); 

E = Eg : (si,a: 2 ,a: 3 ) T ariFi +X2F2 +X3F3. 

( 6 . 6 ) 

Proposition 43. Any element of the Lie algebra of AD + (3) is of the form 
diag(—(ti +<2 + £ 3 ), ti, t 2 , ts), and the corresponding bivector o/A 2 (R 3,3 ) is 

~ 2i^ 1 (® 22 ' + E33') + t 2 (E33' + Eh') + i 3 (Eii' + E22')}- (6.7) 

For example, the isotropic positive dilation diag(l, e 4 , e 4 , e 4 ) is generated by 
—t(En' + E22' + E33'). 

3. Upper triangular shear transformation group: 

An upper triangular shear transformation of £ 3 by vector t = (fi, < 2 , fs) T has 
the following matrix form: 

f 1 

0 

0 

V° 


0 0 0 \ 

1 f.3 t'2 

0 1 h 

0 0 1 ) 


( 6 . 8 ) 


Proposition 44. (6.8) is generated by the following element o/t1 2 (R 3 ’ 3 ): 
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The corresponding element of si (4) is 

/ 0 0 0 0 \ 

0 0 t 3 te-t’p- 

0 0 0 t\ 

\0 0 0 0 ) 


( 6 . 10 ) 


Proof. Let H = a; 3 Ei 2 ' + :ciE 23 ' +a; 2 Ei 3 /. By exp(H) = 1 + H - xix^^yy, 
we get the 6 x 6 matrix form of ^4^* xp( - H ) as following: 


(l 2z 3 

0 1 
0 0 
0 0 
0 0 
\0 0 


2(x 2 + xix 3 ) 
2x\ 

1 

0 

0 

0 


0 

0 

0 

1 

—2x 3 

- 2 ( x 2 - £12:3) 


0 

0 

0 

0 

1 

—2x\ 


o\ 

0 

0 

0 

0 

1; 


( 6 . 11 ) 


The corresponding 4x4 matrix of *SX(4) is 

/l 0 0 0 \ 

0 1 2x 3 2 {x 2 + x 3 x 3 ) 

0 0 1 2xi 

\0 0 0 1 / 


□ 


4. Other triangular shear transformation groups: 

In the ripper triangular shear transformation ( 6 . 8 ), basis vector ei is fixed, e 2 
is sheared towards ei, and e 3 is sheared towards plane eie 2 . The roles of ei, e 2 , e 3 
can be interchanged, resulting in five other branches of shear transformations. 
The six different branches of shear transformations are generated by the following 
six 3-D Lie subalgebras, each being represented by its basis: 


E 12 ', E 13 /, E 2 3 ',’ 

Ei 2 ',Ei 3 /,E 3 2 '; 






E 21 ', E 23 ', E 31 /; 

E 21 ', E 23 ', Ei 3 /; 




( 6 . 12 ) 

E 3 1 ', E 32 ', Ei 2'5 

E 3 1 ', E 32 ', E 2 l'. 







1 

/I 

0 

0 

(A 


For example, the lower triangular shear transformation 

0 

0 

1 

£3 

0 

1 

0 

0 

is 


1 

^0 

£2 

tl 

V 



generated by the following element of yl 2 (R 3,3 ): 

-(tiE 32 ' + t 3 E 2 i' + (O —-^-)E3i'). (6.13) 

5. 3-D Lie subalgebra sl( 2): 

In £ 3 let Xet 4 represent a point, and let 77 £ (R 4 )* represent a plane not 
incident to the point. The subgroup of SX(4)USX - (4) fixing point X and plane 
77 is the group of 2-D projective transformations of the plane with respect to 
the point, and is isomorphic to SL(3) U SL~(3). 
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This embedding of SL( 3) U 777“ (3) into 7>7(4) may be suitable for modeling 
a pin-hole camera. Let C,X £ R 4 , and let 77 £ yl 3 (R 3,3 ) represent a plane, 
then the intersection of line CX with the plane is represented by the bivector 
(CX) • 77 £ yl 2 (R 3,3 ). A null 2-space of 77 represents a unique point on the plane. 
It can be taken as the local representation of the spatial point on the plane. The 
global representation of the point is (CX) A ((CX) • 77) £ A 3 (R 3,3 ). 

In the special case where the point is the origin, and the plane is the one at 
infinity, The 8-D Lie subalgebra s/(3) is spanned by the following elements of 
/1 2 (R 3 ’ 3 ): E(x) A E'(y) for x,y <E R 3 . 

When vectors e 0 ,e 3 are fixed, and the 2-space (ei,e 2 ) is also fixed, the cor¬ 
responding Lie subgroup of 777(3) is SL(2). The Lie subalgebra sl(2) has the 
following bivector basis: E 12 ', E 21 /, En< — E 22 /. 


Lemma 20. The Plucker transform induces the following Lie algebraic isomor¬ 
phism from sl( 4) to so(3,3): 


M = 




e sl( 4) h> 


f —tr(N)I 3x3 + N 
V m o x I 3x3 


where N = (mi m 2 m 3 ), and m,, n 0 £ R 3 for 0 < i < 


—n 0 x I 3x3 \ 
tr(N)l 3 X3 -N T ) ’ 

(6.14) 

3. 


Proof. As M is the derivative of exp(tM) at t = 0, for X, Y £ R 4 , 


(exp(iM)X) A (exp(tM)Y) = (X + 7MX + o(t)) A (X + fMY + o(t)) 

= X A Y + t{ MX A Y + X A MY) + oft) 

= X A Y + t(M 0 I 4x4 + I 4x 4 0 M)(X A Y) + oft), 

(6.15) 

where 

M0l 4x4 (XAY) = iM0l 4x4 (X0Y-Y0X) = ^(MX0Y-MY0X). (6.16) 


So the lift of M £ sl( 4) to so(3, 3) is the 6x6 matrix representation of the action 
of M 0 I 4x4 + I 4x4 0 M upon A 2 (R 4 ) = R 3,3 . For any positive permutation ijk 
of 123, by 


(M 0 I 4x4 + I 4x4 0 M)e 0i = -tr(N)e 0i + moe, + e 0 m„ 
(M 0 I 4x4 + I 4x4 0 M)e jk = -n 0 ei + ejin k - e fe mj, 


we get (6.14). 


□ 
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Definition 19. The following are notations of some linear maps from R 3 to 
3x3 matrices: 

diag : x = a^ei + x 2 e 2 + x 3 e 3 i-g diag(a:i, x 2 , x 3 )\ 


/ 0 x 3 0 \ 

skew : x = X\e 3 + x 2 e 2 + x 3 e 3 i—>• | 0 0 x\ 1 ; 

\x 2 0 0 / (6.17) 

/ 0 0 Z2 \ 

skew T : x = xiei + x 2 e- 2 + x 3 e 3 i-g x 3 0 0 = (skew(x)) T . 

\ 0 X! 0 / 

Proposition 45. The Pliicker transform and the adjoint action of Spin(3, 3) 
induce the following The Lie algebraic isomorphism from sl( 4) to A 2 (R 3 ’ 3 ): let 

N = diag(ni) + skew(n 2 ) + skew T (n 3 ), (6.18) 

and let m 0 , n^ G R 3 for 0 < i < 3, then 

M = ("mo^ N ) € Sl(4) ^ ^{- E ' E '( m o) + EE(n 0 ) - E = E'(m) 

+E<E'(n 2 ) + E>E'(n 3 )}. 

(6.19) 

Proof. The right side of (6.19), denoted by B 2 G rl 2 (R 3,3 ), induces a linear 
transformation X G R 3,3 i—>• B 2 • X G R 3,3 . It has the same matrix form (6.14). 

□ 

Definition 20. Given an n x n symmetric real matrix K, the set of matrices 
M G SL{n) satisfying M T KM = K is denoted by SO (K), called the special 
orthogonal group with respect to K. Its Lie algebra is denoted by so(K). 

Proposition 46. Let K be a 4 x 4 symmetric real matrix of rank > 3. Then 
so(K) is isomorphic to one of so(4), so(3,1), so(3, 0,1), so(2,2), so(2,1,1). 

Proof. Let R G 0(4) be a matrix to diagonalize K by similarity transforma¬ 
tion: R t KR = L, where L is a real diagonal matrix whose nonzero entries are 
the eigenvalues of K. Let T be a real diagonal matrix to change each nonzero 
entry of L to ±1 by matrix congruent transformation: T J LT = N, where N is 
a real diagonal matrix whose nonzero entries are ±1. Up to sign, N is one of 

diag(0,1,1,1), diag(0,-1,1,1), I 4X 4 , diag(-l, 1,1,1), diag(-l,-1,1,1), 

and 

(RT) t K(RT) = N. (6.20) 

Let M G so(K). Then exp(M T )K exp(M) = K. It can be written in the 
following form by using (6.20): 

{(RT)- 1 exp(M)(RT)}' r N{(RT)- 1 exp(M)(RT)} = N. (6.21) 
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So (RT) 1 exp(M)(RT) G SO( N), and so(K) is isomorphic to so(N). 

We start to discuss 6-D Lie subalgebras of type so(K). 

6. so(4): 

It is composed of skew-symmetric matrices, i.e., matrices of the form 

0 x T \ 

-x y x I 3x3 ) ’ 

for x, y G R 3 . The following is a bivector basis of so(4): 


□ 


( 6 . 22 ) 


E 2 ' 3 ' + E 33 , E 3 /p + E 3 i, Ei<2' + E12, 

E T7' XT' B B TT' 

32 ' — ii/ 23 ', H/13/ — H/31', H/ 21' — 

For x = (xi,x 2 ,x 3 ) t and y = (y u y 2 , yz) T of R 3 , let 

Q(x) = a;i(E 3 2' — E 23 ') + £2 (Ei 3 / — E 3 i<) + a; 3 (E 2 i' — Ei 2 '), 
P(y) = yi(E 2 / 3 / + E23) + y 2 (E 3 'i/ + E 33 ) + y 3 (Ei/2' + E12). 


(6.23) 


(6.24) 


If we denote by 



the bivector Q(x) + P(y), then 


xi 

yi 



(xixx 2 + yixy 2 \ 
\ x i x y 2 + yi x x 2 ) 


(6.25) 


where “x so ( 4 )” denotes the Lie bracket of so(4). So so(4) defines its own cross 
product of screw forms. 

If we set = P( x ) + E'(y), then the corresponding trivector of /1 3 (R 3,3 ) 

realizing the cross product of so(4) is 


Ei 2 3 + Ei 2 ' 3 ' + E 2 3 'i' + E 3 p 2 ' = ^"(eo) + -F(ei) + .F(e 2 ) + T(e 3 ). (6.26) 


It corresponds to the matrix diag(l, 1,1,1) of the quadratic form preserved by 
group £0(4). 

7. so(3,1): 

It is composed of matrices of the form 


0 x T \ 
x y x I 3x 3 ) ’ 


(6.27) 


for x,y G R 3 . It has the following bivector basis in H 3 (R 3,3 ): 


E 3 ' 3 ' 

— E 23 , 

E 3 u' 

— R31 5 

E]/ 2 ' — E 12 

E 3 2 ' 

— E 23 ', 

Ei 3 / - 

- E 31 /, 

E 2 I' ~ E 12 ' 


(6.28) 


Let 

Q(x) = (E'E' - EE)(x), P(y) = (E>E' - E<E')(y). (6.29) 
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If we denote by 



the bivector Q(x) + P(y), then 


X 1 

yi 



/xixx 2 -yixy 2 \ 
V x i x y 2 + yi x x 2 y 


(6.30) 


where “x so ( 3> i)” denotes the Lie bracket of so(3,1). So so(3,1) defines another 
cross product of screw forms. 

If we set = E(x) + E'(y), then the corresponding trivector realizing 

the cross product of so(3,1) is 


—Ei 23 + Ei 2 / 3 / + E 23 /i/ + E 3 i/ 2 / — —J-'(eo) + .F(ei) + J r (e 2 ) + .F(e 3 ). (6.31) 


It corresponds to the matrix diag(—1,1,1,1) of the quadratic form preserved by 
group SO( 3,1). 

8. so(3,0,1): 

It is composed of matrices of the form 


0 x T \ 

0 y x I 3x3 J ’ 

for x, y £ R 3 . It has the following bivector basis in yl 3 (R 3,3 ): 


Let 


E 23 , 

E 32 / — E 23 /, 


E 3 i, 

Ei 3 / — E 3 u, 


Ei 2 , 

E 2 i/ — Ei 2 /. 


If we denote by 



= EE(x), P(y) = (E>E , -E<E , )(y). 
the bivector Q(x) + P(y), then 


(6.32) 


(6.33) 


(6.34) 


Xl 1 v ( x 2 

yi X^Ol) y 2 


xi x x 2 

X 1 X y 2 + yi x x 2 


(6.35) 


where “x so ( 3i0 ,i)” denotes the Lie bracket of so(3,0,1). So so(3,0,1) defines the 
same cross product of screw forms with that of se(3), and the two algebras are 
isomorphic under the transpose of 4 x 4 matrices. 

The trivector realizing the cross product of so(3,0,1) is 


Ei 2 ' 3 ' + E 23 'i< + E 3 i' 2 / — -E(ei) + J r (e 2 ) + .F(e 3 ). (6.36) 

It corresponds to the matrix diag(0,1,1,1) of the quadratic form preserved by 
group SO( 3, 0,1). The transpose of matrices from se(3) to so(3,0,1) is realized 
in A 2 (R 3 ’ 3 ) by A 2 J. 

9. so(2,2): 



56 


Hongbo Li, Lei Huang, Changpeng Shao, Lei Dong 


It is composed of matrices of the form 


/ 0 

X\ 

X 2 

X 3 \ 

-Xi 

0 

2/3 

-yi 

x 2 

2/3 

0 

2/1 

\ X 3 

-yi 

—2/1 

0 / 


(6.37) 


It has the following bivector basis in yl 3 (R' 


3,3V 


E2'3' + E 23 , 
E 32 ' — E 23 ', 


E 3 'i' — E 31 , 

E 13 / + E 31 /, 


El'2' — E12, 

E 2 I' + E 12 '. 


Definition 21. For x = (x\, x 2 , x 3 ) T and y = (j/i, yi, ys) T of] 


X 2 Vz - x 3 y 2 \ 


'-(x 2 y 3 - x 3 y 2 )\ 


x 1 y := ~{x 3 yi~xiy 3 ) , x x 2 y := -(x 3 yi + X\y 3 ) 


^{xiy 2 - x 2 y\) J 


xiy 2 + x 2 yi ) 


In so(2,2), for x = (xi, x 2 , x 3 ) T and y = ( 2 / 1 , y 2 , 2/ 3 ) T ofR 3 , let 

Q(x) = Xi( — E32' + E23') + X 2 ( E13/ — E31') + X3(E2i' + E12'), 
P(y) = J/l(E2'3' + E23) + J/2(E 3 '1' — E31) + 2/ 3 (Ei'2' — E12). 


(6.38) 


(6.39) 


(6.40) 


If we denote by 


x 

y 

Xl 

yi 


the bivector Q(x) + P(y), then 

x so(2,2) 


xi xi x 2 +yi x y 2 
yi x 2 x 2 - y 2 x 2 xi 


(6.41) 


where “x so (2,2)” denotes the Lie bracket of so(2,2). So so(2,2) defines its own 
cross product of screw forms, 
x 

,y. 


If we set 


= E( x) + E'(y), then the corresponding trivector realizing 
the cross product of so( 2 , 2 ) is 

E 123 — Ei 2 ' 3 ' + E 23 U' + E 3 ]/ 2 ' = — -F(eo) — ^"(ei) + ^(e 2 ) + -F(e 3 ). (6.42) 

It corresponds to the matrix diag(—1, —1,1,1) of the quadratic form preserved 
by group 50(2, 2 ). 

10 . so( 2 , 1 , 1 ): 

It is composed of matrices of the form 


(6.43) 


( 0 

Xl 

X2 

x 3 \ 

0 

0 

2/3 

-yi 

0 

2/3 

0 

2/1 


yi 

—2/1 

0 / 


It has the following bivector basis in /1 3 (R 


3,3y 


E23, 

E32' — E23' 


—E 3 i, 

Ei 3 / + E 3 i/ 


E12, 

E 21 ' 


E 


12 ' 


(6.44) 
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In so(2, 1,1), for x = (x 1 ,x 2 ,x 3 ) T and y = {yi,y 2 ,y 3 ) T of R 3 , let 


Q(x) — Xi (E 32 ' — E 23 ') + x 2 (E 13 / + E 31 /) + a; 3 (E 2 i' + E 12 '), 

P(y) = 2/1E23 — y 2 E 3 i + 2/3E12. 


If we denote by 



the bivector Q(x) + P(y), then 


Xl 1 v X2 

yi ) x so(2, 1 ,i) ( y 2 


xi Xi x 2 

yi x 2 x 2 - y2 x 2 xi 


(6.46) 


where “x so ( 2 ,i,i)” denotes the Lie bracket of so(2,1,1). So so(2,1,1) defines its 
own cross product of screw forms. 

If we set = E(x) + E'(y), then the corresponding trivector realizing 

the cross product of so( 2 , 1 , 1 ) is 


Ei2 / 3 / + E 23 U' + E3i'2' — —.F(ei) + J 7 ^) + ^(es). (6-47) 

It corresponds to the matrix diag(0, —1,1,1) of the quadratic form preserved by 
group SO( 2,1,1). 

Definition 22. The Hadamard product of vectors in R 3 is the following multi¬ 
linear, associative and commutative product: 


( Xl\ f X2 \ / XlX 2 \ 

yi © I 2/2 : = yiy? ■ ( 6 - 48 ) 

) \ z 2 ) \ Z\Z 2 ) 

The i-th Hadamard power of a vector x £ R 3 is 

©*x:=x©‘--©x (i times). (6.49) 


The exponential function of x in the Hadamard product is denoted by e 0x . In 
particular, 


©°x = e 00 



(6.50) 


11. 6 -D Lie subalgebra of general anisotropic dilation group: 

The general anisotropic dilation group is composed of matrices of the form 


(1 

0 

0 

°\ 

tl 

Ai 

0 

0 

£2 

0 

A 2 

0 


0 

0 

A 3 / 


The matrix represents either an anisotropic dilation centering at an affine point, 
or a pure translation. The corresponding Lie subalgebra is denoted by gad( 3), 
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and is composed of matrices of the form 


/ 0 0 

Ml Hi 

m 2 0 
\u 3 0 


0 0 ^ 
0 0 
H2 0 
0 Ms/ 


(6.52) 


Its bivector basis is spanned by elements of the form E'E'(x), E=E , (y) for 
x, y e R 3 . 

Let 

Q(x) = xi'Eyy + Z 2 E 3 /P + a; 3 E 1 / 2 ', 

P(y) = 2/1F1 + 2/2F2 + 2/3F3, 


(6.53) 


where the F, ’s are defined by (6.5). If we denote by ^ ) the bivector Q(x) + 
P(y), then 


£ M* 


xi © y 2 - yi © x 2 

0 


(6.54) 


where u x ga d” denotes the Lie bracket of this Lie algebra. 

12. 7-D invariant group of R 2,2 and its 6 -D subgroups: 

Let L 4 be a fixed 4-space of signature R 2,2 , and let its orthogonal complement 
in R 3,3 be L 2 . Then L' 2 has signature R 1,1 , so its null 1-spaces represent a pair of 
non-intersecting lines in space, say li and 1 2 . The null 1-spaces of L 4 is denoted 
by 7V(L 4 ); it represents all the lines in space incident to both li and 1 2 . 


Proposition 47. Any projective point on line li or 1 2 is on infinitely many lines 
of N( L 4 ), while any other projective point is on one and only one line of NfiLfi). 
Similarly, any plane passing through li or 1 2 also passes through infinitely many 
lines of N( L 4 ), while any other plane passes through one and only one line of 
N( L 4 ). 

Proof. Let X 3 be a null 3-space of H(R 3,3 ). When X 3 represents a point of 
li, then it is incident to all the lines connecting X3 and 1 2 ; when X3 represents 
a plane through li, then it meets line 1 2 at a point Y3, and so is incident to all 
the lines connecting Y 3 and li. 

When X3 represents a point on neither li nor 1 2 , then in R 3,3 , the intersection 
of the 3-space X 3 with the 4-space L 4 is a null 1-space or 2-space. If point X 3 
is on more than one line of IV(L 4 ), then li,l 2 must be coplanar, violating the 
requirement that the two lines do not intersect. So point X3 is on on and only 
one line of !V(L 4 ). Similarly, when X3 represents a plane supporting neither li 

nor 1 2 , it must be incident to one and only one line of 7V(L 4 ). EH 

Consider the subgroup of 50(3,3) that leaves the 4-space L 4 invariant. De¬ 
note the subgroup by /nu(L 4 ). Then 


Inv{ L 4 ) = (50(1,1) © 50(2,2)) U (50"(1,1) © 50"(2,2)). (6.55) 


It is a 7D Lie subgroup. 
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A typical example is that li is an affine line, say li = eoi, the line through 
the origin and in direction ei, and I 2 is a line at infinity, say I 2 = e 23 , the line 
at infinity normal to ei. Then 7 V(L 4 ) is composed of all lines meeting li and 
perpendicular to it. Then L 4 = E 22 ' 33 ', and denote 


Inv(2,2) := Inv(E 2 2'33')- (6.56) 

The Lie algebra inv( 2,2) of Inv( 2,2) is spanned by Eir, E 2 2 ', E 33 /, E 23 , E 2 ' 3 ', 
E 23 ', E 32 '. Its matrix form in si (4) is 


A 0 
0 D 


(6.57) 


where A, D are two 2x2 matrices such that tr(A) + tr(D) = 0. The group 
Inv( 2 , 2 ) is composed of matrices of the form 


A 0 
0 D 


or 


0 B 

C 0 


(6.58) 


where det(AD) = 1, and det(BC) = —1. 

For a general 4-space L 4 of signature R 2 ’ 2 , Inv(Tn) is isomorphic to (GL(2) x 
GL( 2))/*, where two elements are equivalent if and only if they differ by scale. 
The space L 4 is called the space of coupled projective screws with respect to the 
pair of axes li, I 2 . The projective screw ratio refers to det(A)/ det(D) in the first 
case of (6.58), and det(B)/ det(C) in the second case. 

A typical 6 -D Lie subalgebra of inv( 2, 2) is spanned by E 2 2 ', E 33 /, E 23 , E 2 ' 3 ', 
E 23 ',E 32 '. For x = (x 1 ,X 2 ,x 3 ) t and y = (yi, y 2 , y 3 ) T , let 

P(x) = a;i(E 33 / -(- E 2 2 ') + a" 2 E 2 3 + x 3 E 2 ' 3 ', , ^ 

Q(y) = Vi (E 33 / — E 22 ') + 2 / 2 E 32 ' + y 3 E 2 3 '. 

Then 

P(x) X Q(y) = 0, 

P(x) x P(y) = -{x 2 y 3 - ^’ 32 / 2 )(E 33 / + E 22 O + 2(xiy 2 - x 2 yi)E 23 

+2(x 3 y± - xiy 3 )E 2 ' 3 ', (6.60) 

Q(x) x Q(y) = {x 2 y 3 - x 3 y 2 )(E 33 > - E 22 ') + 2(xiy 2 - x 2 yi)^ 3 2 ' 

+ 2 (a; 3 7 /i - xi 2 / 3 )E 23 '. 


For the whole Lie algebra sZ(4), by (6.17), any element of it is naturally 
decomposed into the direct sum of five 3-D vectors: 

—tr(diag(c)) n T 

m diag(c) + skew(u) + skew T (d) 

where m, n, c, u, d £ R 3 . They allow for defining a matrix of sZ(4) as a “super¬ 
screw” in R 3 x R 3 x R 3 x R 3 x R 3 , and defining the Lie bracket as a “supercross 
product”. Such a representation of si (4) is called the superscrew representation. 
For symmetry consideration, instead of choosing the 3-D Lie subalgebra decom¬ 
position (6.61), we make the following decomposition of sl( 4) into 3-spaces: 


(6.61) 
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Definition 23. When sZ(4) is represented by bivectors of H 3 (R 3,3 ), a superscrew 
of 3-D projective geometry is defined as follows: for any x t £ R 3 where 1 < i < 5 , 

EE(xi) + E'E'(x 2 ) + E = EJ,(x 3 ) (6.62) 

+(E>E' + E<E')(x 4 ) + (E>E' - E<E')(x 5 ). 

Definition 24. The symmetric cross product of two vectors of R 3 is defined as 
follows: 

( X 1 \ / x 2 \ / y\Z 2 + Ziy 2 \ 

I Vi I * I V 2 I = I z 1 x 2 +x 1 z 2 I . (6.63) 

\zij \ z 2 J \a;ij/ 2 + yix 2 ) 

By direct computation, we get 

Proposition 48. The Lie bracket of si(4), when represented in the superscrew 
form, is the following: for Xj,yj £ R 3 , 




where 


zi = -xi © y 3 + x 3 0 yi +xi *y 4 - x 4 *yi +xi x y 5 + x 5 x y 4 , 

Z 2 = -x 2 * y 4 + x 4 * y 2 + x 2 © y 3 - x 3 © y 2 + x 2 x y 5 + x 5 x y 2 , 

Z 3 = (~xi • y 2 + x 2 • yi)e 00 - x 4 © y 2 + x 2 © y 4 + 4(x 4 © y 5 - x 5 © y 4 ) x e 00 , 

z 4 = J(x 4 *y 2 - x 2 *y 4 ) - (x 3 x e 00 ) ©y 5 +x 5 0 (y 3 x e 00 ) 

x 4 x y 5 - x 5 x y 4 , 

z 5 = -^(xi x y 2 + x 2 x y 4 ) + (x 3 x e 00 ) © y 4 - x 4 © (y 3 x e 00 ) 

+X 4 x y 4 + x 5 x y 5 . 

(6.65) 


7 Conclusion 

In this paper, we establish a rigorous mathematical foundation for the line ge¬ 
ometric model of 3-D projective geometry. We also extend screw theory from 
rigid-body motions to projective transformations, in hope that non-Newtonian 
mechanics may find an algebraic language for the development of virtual work 
of projective motions. 
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The connection between the Cl( 3,3) model of 3-D projective geometry and 
the Cl{ 4,1) model of 3-D conformal geometry is an interesting topic, and will be 
investigated in another paper. This connection makes it possible to investigate 
3-D non-Euclidean geometry together with its various realizations in Euclidean 
space via pin-hole cameras of R 4,1 . 
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